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CHAPTER I 


INTRODUCTION 


An important development in recent collision experiments has 
been the evidence of autoionizing compound negative ion states with 
lifetimes from ome to vomual Seconds; considerably shorter than the 
lifetimes of electronically excited states (~10 ~ sec.). These 
compound states have been used to explain resonances in elastic 
scattering cross sections (Simpson and Fano, 1963), the onset of 
inelastic processes where no states of the target exist (Schulz end 
Philbrick, 1964; Chamberlain, 1965), forbidden transitions (Baranger 
and Guerjoy, 1957), and dissociative attachment (Schulz, 1958). 

Most of the states observed in electron scattering experiments 
have energies lying slightly (0.1 - 1.0 eV) below known excited states 
of the target and the usual decay process is autoionization. Some 
Specific examples of targets are hydrogen, helium, and the hydrogen 
molecule. 

In hydrogen, the first excited state is “s 26°09 <p ZO aie 
10.204 eV above the ac ls ground state. A resonance has been observed 
by Schulz (1964) at 9.75 4 .15 eV, just below the 2s or 2p energy. 


At’ 19.3 eV above the ground state of the helium atom there is a 


Sharp He™ resonance which is so well known that it is used for energy 











eolemealibration (senulzs 1964; Simpson, Kiyatt, et als, 19654: and 


—S hols) staves C1 meltum ay 


many others). This is just below the 
an everage energy of 20.2 eV. Fano and Cooper (1965) have also dis- 
cussed higher He™ resonances, such as the <p (2s)"2p A pt ley “anid 
> 2s(sp)* at 58.2 eV. These lie just below the tg (2s)* and 

3p esep states at 5/.9 and 58.4 eV respectively. 

The pre-ionization states of helium can be considered as reso- 
nances of Het plus an electron. Madden and Codling (1964) have 
observed a set of these autoionizing levels, and Cooper, Fano and 
Prats (1963) have shown that they correspond to a series (2snpins2p) 'P, 
fie 2,5, -... which converges Lo Het 2p (or 2S) Bile, HOS oes 

Kuyatt, Simpson and Mielczarek (19650, 1966) have observed pairs 
of reSonances starting about .9 eV below the PA, states of the 
hydrogen molecule, and studied the isotope effect in HD and D,,- In 
these experiments, spacings of levels were obsServed which correspond 


closely to the vibrational levels of the excited H, (HD,D,) states, 


9) 
which seems to indicate that the third electron in the compound state 
became trapped in the field of the excited target state, but far 
enough away that it had little effect on the nuclear motion. Golden 
ae cel (1965) report substantially the same result. 

Various methods have been proposed and used for calculating the 
erergies and wave functions of these compound states, particularly 
for scattering of electrons off atoms. The problem is that the normal 


veriation method is inapplicable, as an infinity of states of the 


Same symmetry and equal or lower energy exist for each reSonance 





Solsisuane, Of the ground state or the target plus a iree electron with 
the proper angular momentum and a continuum of energy values. 

Fano (1961) and Burke and Schey (1962) have developed methods 
which actually deal with the scattering cross section. The Fano 
theory is particularly useful for analysis of line shape, while Burke 
and Schey use a close-coupling technique to calculate the cross section 
and study resonances as part of the total scattering problem. 

Holgien (1958a,b), Holgien and Midtal (1965, 1966) Mandl and 
Herzenberg (1963, series), O'Malley and Geltman (1965), Taylor and 
this author (1965), and Lipsky and Russek (1966) have devised quasi- 
stationary methods for calculation of the energy and wave function of 
resonant states without predicting any of the features of the cross 
section except the approximate position of the resonance. Of the 
above, Taylor and Williams (1965, hereinafter denoted as TW) were the 
first to do calculations on systems other than two-electron atoms. 


The Taylor-Williams calculation was on the diatomic molecule H, plus 


2 
an electron, and the results compared quite well with the experiments 
of Kuyatt, et al. (1965) and Golden and Bandel (1965). Since that 
time, Mandl and Herzenberg (1965) have done a parametrized calculation 


on the ie state of H. plus an electron with their modified Kapur- 


C 

Peierls technique. 
Taylor, Nazaroff, and Golebiewski (1966, hereinafter denoted as 

TNG), with the aid of the author, have presented a formal justification 


of the physical quasi-variation method of TW, and presented a com- 


plete quasi-stationary method for calculating the energy of resonances 





which includes the TW method as a Special case. This work also lists 
and classifies three types of resonances, including the examples 
above and two other types which will be discussed later. 

Therefore, one of the aims of this work is to present the 
Taylor-Williams quasivariation method with its original physical 
ereuments, and the formal justification of it as presented by TNC. 
This is done in Chepter II. 

The second aim is to apply these methods to small diatomic 
molecules (in particular, H,); calculating resonant states and com- 
paring to electron scattering experiments (Chapter III, ff.). 

Pinel, im-econjunetron with Lhe Second aim, the dravomie 
molecule calculation methods of Toe (1960) Taylor and Harris (1963, 
series), and Taylor (1964) ‘or bound states of six electrons or less 
have been revised to include the quasi-stationary methods, and have 
been extended to allow both bound and quasi-stationary calculations 
of energies and wave functions of larger diatomic systems. The 
latest revision, presented and discussed in the Appendices, can do 


calculations on up to twenty electrons. 





CHAP Ee 


THEORY 


The formation, lifetime, and decay of the compound states 
mentioned above can only be described rigorously by wave packets 
composed of combinations of bound and continuum states of the Ntl 
electron system. The reSonance is Simply a highly localized wave 


Le to ey seconds, so that with the 


packet which remains so for 10. 
proper reStrictions, it can be treated as stationary during this time, 
This is similar to the treatment of electronically excited states as 
Stationary by ignoring the light field interactions, where inclusion 
of the light field leads to decay in approximately 10° seconds. 
Without any restrictions, a carefully executed variational cal- 
culation will lead to the conclusion that no compound states exist, 
Since one can construct an infinite number of states of the same 
symmetry and equal or lower energy. Thus, even to apply the Hylleraas- 
Undheim procedure, treating the resonance as the aga excited state, 
one would have to diagonalize an infinite secular determinant, taking 
the uppermost root as the resonance energy. 


The original work of Taylor and Williams (TW, 1965) presented a 


quasi-variation method based on physical arguments for selecting the 





mROoven MesurictLions so that valad results could still be obtained. 

Briefly the method stated that if an N+l electron bound state 
function could be guessed such that it nearly represented the physical 
situation, it would also closely resemble the proper linear com- 
bination of wave packets, and the matrix elements of interaction of 
ieee DOUG function with the infinite number of continuum functions ~‘* 
would be so small compared to the bound-bound interactions that the 
infinite Hamiltonian matrix would be nearly numerically block diago- 
nalized, separating the bound state in the continuum from its adjacent 
continuum. 

The rules given for selecting these Ntl electron functions 
were that they must be doubly excited, autoionizing configurations, 
EMtewith small autoionization probability in order to obtain Long 
lifetimes. The small autoionization probability could be assured by 
selecting "weakly correlated" orbitals (i.e., orbitals having small 
1/r, matrix elements) for the trial wave function. This choice of 
Tunetions causes both the Hamiltonian and overlap matrices to be 
nearly block diagonal in the prescribed manner, since the double 
excitation produces an excited core which is orthogonal to the ground 
state of the target. Integration over the inner electrons will make 
the bound-continuum matrix elements small. Hence the total secular 
determinant is approximately block diagonalized into a finite part 
mixing multiply excited configurations of the compound state and an 
infinite part which can be ignored. 


These arguments work quite well for certain types of resonances, 





but considerable experimental evidence has appeared that there are 
other resonances which the quasi-variation method would not be able to 
calculate. Taylor, Nazaroff, and Golebiewski (ING, 1966) analyzed 

the experimental data and proposed that it could be explained by 

three types of resonance, which they classified as Type I Core-Excited, 
Type II Core-Excited, and Single Particle (hereinafter denoted by 

CEl1, CE°®, and SP, respectively). 


The resonances of He, H, and H, discussed above, which lie 


PB: 
O.1 to 1.0 eV below excited states of the target, are CEl. They could 
be formed by the virtual process of the impinging electron exciting 
one of the electrons of the target and becoming trapped in the slight 
potential well left when the excited electron moves out of its way, 
with a resultant binding energy of 0.1 - 1.0 eV. The CEl resonances 
are the ones which can be handled by the quasi-variation method, and 
the excitation of one of the target electrons explains the classi- 
fication Core-Excited. 

The potential well can be considered as being formed by 
additional amount of nuclear charge that the impinging electron sees 
when the excited core electron moves away. 

One of the examples in many textbooks (for instance Messiah, 
1061, Vol. I, p. 88) is the case of virtual levels of the square well 
potential, in which it is shown that if the walls of the finite 
CF = ©) well were made finite, “long-lived" localized wave packets 
(bound states in the continuum) can be formed only at those energies 


in the continuum of the well (E > V_) at which bound states would have 





existed had ve not been reduced to a finite value. 

Based on this argument, TNG proposed the existence of C2 
resonances as virtual states of the potential well caused by core 
excitation. These resonances should be quite broad (short lived - 
width is inversely proportional to lifetime) compared to CEl, since 
they are adjacent to the continua of both the ground state of the 
target and the excited state of the target, and have no binding 
energy to provide stability. Experimental examples cited as indicative 


4 


that CE2 resonances exist were the Hg™ (Sdeosoo.) Pp seen by 


By 2 
Kuyatt, Simpson, and Mielczarek (1965) at 4.89 eV, and an elusive 
peak seen by the same experimenters in the elastic scattering of 
electrons from helium at 21.4 to 21.6 eV; JUSt above tne n= 2 levels. 
Schulz and Philbrick (1965) and Chamberlain and Heideman (1965) 
observed a stronger peak in the inelastic cross section of He corre- 
sponding to nearly the same energy. TNG attribute this resonance to a 
CE2 above the n = @ level of He, rather than a CEl below the n = 3 
level, as the necessary binding energy below the n = 3 level would 
indicate a longer lived resonance and therefore a narrower peak than 
that observed, and an electron affinity of 2 eV for the n = 3 level as 
epposed £0 .1 to .6foer then = 2 level, 

In the same spirit, it is argued that if an excited target state 
can support virtual levels, it should be possible for the ground 
state to do so also. Several experimental observations exist which 


cannot be explained except by postulating such states, which TNG 


classified as Single Particle since the physical model assumed is the 





nearly unperturbed ground state plus an orbiting planetary electron. 
Miese states are also expected to be quite broad. 

The experiments leading to these conclusions are the N,, - 
electron scattering experiments of Schulz (1962, 1964, 1966), the aes 
electron dissociative attachment experiment of Schulz and Asundi (1966), 
end structure in the elastic He - electron scattering cross section 
at about 0.45 eV observed by Schulz (1965). 

The introduction of the CEH2 and SP resonances requires a more 
effective procedure for calculating the energy of resonances. 

TNG, with the assistance of the author, presented a stabilization 
method which generalizes the quasi-stationary methods mentioned above 
(page 3), and incorporates a way of handling CKe and SP as well as 
CEl resonances. 

The methods of calculating the resonant energies are vasically 
refinement methods. The wave function of the resonant state is 
written as a superposition of configurations selected on the basis of 
experimental results and intuition. The refinement method is used to 
improve the wave function by calculating correction terms. Even 
though this method is an approximate method, it will be justified 
below on formal grounds. The most direct justification of the method, 
however, is the fact that it allows one to calculate ab-initio very 
good approximations to the experimentally observed resonant energies. 

In order to gain perspective it is well to consider for the 
moment how one might calculate an approximate wave function and energy 


of an excited bound state of an atom or molecule. If the State in 





10 


question is the Mth excited state of a given symmetry one might choose 
a large number of properly symmetrized basis functions (greater than 
M), diagonalize the Hamiltonian of the system and optimize the Mth root 
according to the Hylleraas-Undheim procedure. However, if the state is 
very highly excited such a procedure will in general not give a good 
meowlt. One is thus forced to try other methods. There are two 
possible alternative approaches: a variation-perturbation method or a 
quasi-variation technique. In a variation-perturbation method one 
@rooses a zero-order wave function and then corrects it by variational 
versions of perturbation techniques. Such perturbation methods include 
the inverse nuclear charge method (Knight and Scherr, 1963), the method 
of Kirtmann and Decious (1966) and the methods of Hirschfelder et al. 
(1965) and Byers Brown et al. (1966). These methods are convenient to 
oeemeemerne System 1S nol too large. For larger systems One 15 

forced to uSe quaSi-variational approaches in which a limited number 
@pecontiourations is Selected on intuitive grounds. Dhese coniig- 
urations are then optimized variationally by means of linear and non- 
linear parameters. In such methods one must be careful that the final 
wave function does have dominant those configurations which one feels 
ene ineuitive grounds should be dominant. Ome adds other confisurations 
to the wave function to See if the mixing 1S appreciable and constartly 
checks whether the energy is stabilized. In other words, in practice 
one must uSe intuition to a Large degree to guide the variational 
method towards the deSired result. 


All of the methods mentioned above depend for their success on 





it 


the choice of trial functions. The Hylleraas-Undheim method will give 
poor results if the basis set is chosen poorly; the perturbation 
methods will give poor or entirely wrong results if the zero-order 
function is not chosen properly; the quasi-variation of linear and 
non-linear parameters will also give poor or incorrect results if the 
Pemtieuralions chosen do not fit the physics. In the following 
discussion it will be shown that the situation is much the same in 
resonent-state calculations. 

In order to discuss resonances let us consider a scattering 
experiment in which a projectile is scattered by a target. The target 
might be a nucleus and the projectile a neutron or proton as in nuclear 
figstcs, Or the target might ce an atom or flea and the projectile 
om clectron as in electron-impact spectroscory. Let us consider an 
eeperiimemt in which a beam of projectiles is fired from @ pun ab the 
target, scattered by the target and then detected in a detecter. The 
entrance channel is thus an incoming projectile with the target in the 
ground state. If the energy of the projectile is insufficient to 
excite the target, the scattering is elastic and the exit channel will 
be an outgoing projectile with the target in the ground state. Far 
from the scattering center (i.e. in the detector) the wave function 
‘or the total system is thus a product of the ground-state target wave 
function and a sum of incoming projectile waves and outgoing spherical 
projectile waves. 

A traditional way of describing an elastic resonance in nuclear 


physics (Blatt and Weisskopf, 1952) is to say that at certain kinetic 





Bas 


ener;zies the incoming projectile peretrates the target and forms a 
eeeounid syste:.. The compound system exists for a short time and then 
decays vack into a tarret in the original state and an elastically 
Scattered electron. The presence of the compound state is evidenced 
Oe sharp and narrow resonance in the elastic scattering cross 

section at the so-called resonant energy. The possible resonant 
energies, i.e. the energies of the possible resonant states, are ob- 
tained by excludire the incoming and outgoing waves and solving for the 
energy levels ot the compound projectile-tarret system as if it were 

a oound system. 

One possible way of excluding the incoming and outgoing continuum 
contributions to the resonant wave function is to put the compound 
system in a box. In nuclear physics this is a reasonable procedure. 
The surface of the box would be the surface of the compound nucleus. 
Since nuclear forces are short-ranged, one can meaningfully define a 
Mielear Surface and can partition physical space anto™imside the 
compound nucleus and outside the compound nucleus. However, in the 
ease of electron-atom or electron-moiecule scattering it is difficult 
to define a vox because of the long-range Coulombic electron-electron 
forces. One achieves the effect of putting the compound system in a 
vox by partitioning the Hilbert space of the compound-system Hamiltonian 
instead of physical space. The Hilbert space is partitioned into 
P-space which contains ali functions sexpandable insterms of the 
asymptotic exit channels and into Q-space which contains all functions 


having no projection on the asymptotic exit channels. Hence, one 





13 


diagonalizes the compound-state Hamiltonian in Q-space. The solutions 
so obtained can have meaning only in the situation in which the 
electron is confined to the vicinity of the target. The electron has 
thus been put in a box whose radius is the range of the effective 
scattering potential without even specifying that radius. It will 

be shown formally below that diagonalizing the Hamiltonian in Q-space 


is equivalent to diagonalizing the pseudo-Hamiltonian QHQ = in the 


"QQ 
full Hilbert space. The operator Q is a projection operator which 
selects Q-space. The method of calculating resonances thus reduces 

to diagonalizing the compound-system projectile-target Hamiltonian 

with as large a basis set as is possible which does not include the 
asymptotic exit channels. 

Let us consider the case of electron-atom scattering. The 
projectile is an electron and the target is an N-electron atom presumed 
to be in its ground state. The total energy of the system is the 
kinetic energy of the incoming electron and the ground-state energy of 
the atom. In order to obtain the resonant states (and energies) of 
the (N + 1)-electron atom one looks for apparent bound states of the 
(N + 1)-electron Hamiltonian in the vicinity of the experimental energy. 
Therefore, in order to calculate the resonant energies of the electron- 
atom system, one must diagonalize the full (N + 1)-electron Hamiltonian 
ot the total system as mentioned above. 

The problem which is encountered is that as mentioned above the 
energy of the reasonant state will not in general correspond to the 


lowest root of the (N + 1)-electron Hamiltonian. The situation is 





a 


thus similar to that encountered in solving for the energy of a very 
highly excited bound state, hence one can use either perturbation- 
variation or quasi-variation methods to zero in on and stabilize the 
@eeired root. In order for a root-stabilization method to work the 
energy of the intuitively-guessed wave function must be close to the 
actual resonant energy. It is only in such a case that one can safely 
discard any configuration which seriously alters the energy. One 
expects a root-stabilization method to work because the electrons most 
involved in the resonances cannot be highly correlated. If they were 
highly correlated, the resonant state would not possess a long enough 
lifetime to be observed experimentally. The low correlation indicates 
iaoea one-electron picture (i.e., trial function as product of one 
electron orbitals, or small linear combination of products) is 
meaningful. One therefore expects the energy of the resonant state 

to have essentially settled down to its final value, i.e. to have 
become Stabilized, after mixing a few configurations with the intui- 
tively guessed one(s). 

A word of caution is necessary on the uSe of a stabilization 
method. If a root cannot be stabilized with a small number of configu- 
rations either an incorrect wave function has been chosen to start 
with or the basis set is not a good one. If very closely-spaced roots 
Derin to occur, then the concept of root-stabilization becomes 
meaningless and the method is inapplicable. This would be the case if 
the resonant state belonged to a group cf very closely-spaced resonant 


States. The desired state would be a member of a group of almost 





i 


degenerate resonant states in this case. In such a case one would have 
to consider all of the almost degenerate resonant states as a unit. 

The quasi-variation method of Taylor and Williams (1965) and the 
equivalent method of O'Malley and Geltman (1965) can now be seen to 
be incomplete. In these methods the Q-space was constructed by 
eliminating from the basis set all configurations containing those 
target states included in the exit channels, i.e. they eliminated all 
parent target states lying below the resonant state. Therefore, these 
methods could not be used to calculate single-particle or type II 
core-excited resonant states because they lie above their parent 
target states. These methods do give results for core-excited type I 
resonant states comparable with experiment eee the core-excited 
resonant states do not involve to a large degree the lower lying parent 
target states. Because of the small correlation between the extra 
electron and the target, the ground target states do not mix very 
strongly. In other words, the quasi-variation method is a way of 
zeroing in on the root for those resonant states in which the wave 
function does not include appreciable overlap with the initial target 
states (i.e. those target states occurring in the exit channels). The 
Q-space of the quasi-variation method is a subspace of the more 
reneral Q-space of the stabilization method. 

The obvious fault of the stabilization methods is that they do 
not calculate lifetimes. One could stabilize, especially in the case 
of a single-particle or type II core-excited resonant state, on a root 


corresponding to a resonant state which is too short-lived to be 





16 


riserved experimentally. The perturbation method will definitely 
merculate stuct. short-lived resonant states if the choice of zero-order 
feve function forces it to do so. In the case of type I core-excited 


meee art, states we are confident that the lifetimes are of the order 


13 


uM 


ong te OnE Om seconds. For single-particle of type II core- 


excited resonant states experience and intuition are not well enough 
develoned to <zllow even rough estimates. 

Now that the quasi-stationary methods have been described, it 
remains to show that the eigenfunctions and eigenvalues they yield do 
correspond to resonant states. In this section the Feshbach theory 
of resorant reactions (1956, 1962) will be used to establish this 


ccrrespondence by showing that eigenfunctions of QHQ = do lead 


“aq 
eee breiuv-Wirser form in the cross section and correspond to a4 
decaying bound state. 


The equation to be solved is the (N + 1)-electron Schrédinger 


equation for the projectile-target system 


HY(2y 29 ++ Tua) = EY (2525) +++92y42) (1) 


for & given energy E. The operator H is the full Hamiltonian of the 


system and ¥(2 


PyoEoo+ + oZy47) 1s the total wave function. The procedure 


used by Feshbach (1962) is to partition space with two projection 
operators P and Q such that they satisfy the relationships 


PY 2, 225) oe Zy41) eso ¥(2y To? eee oe) (2) 
eal Gell 





IE 


and 
Par7Gr= i 
ee a ug (3) 
PQ = QP = 0 


The projection operator P projects out of the wave function all com- 
ponents on the available asymptotic exit channels at the given energy 
E. 

The Schrodinger equation, Eq. (1) can now be solved formally 
in a way which will show the resonances explicitly. Writing 


Y = QY + PY Feshbach partitioned Eq. (1) into a set of two coupled 


equations 

(Hpp-E ) BES =O (ha) 
and 

(Hog-E )a¥ = -HopPt (4b) 


where H,, = PHP, etc. Solving Eq. (4b) for Q¥Y and inserting the 


result into Eq. (4a) yields 


1 
(ee ae "op PY = EPY (Sa) 
QQ 
or 


(App + ee = EPY (5b) 


where Toate is the optical potential. In order explicitly to introduce 


the resonances into the picture one must consider the equation 


(Faq : c,,), =o . (6a) 
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If this equation has eigensolutions, they must be eigenfunctions of 
Q with eigenvalue unity, i.e., Qs. = o? and Q@ can be expanded in 


terms of its eigenfunctions as 
Tey (6b) 
n 


where the summation includes the integral over the continuum if 


necessary. Inserting this expression into Vo yields 


pt 
Hpg|8n)(8n| Hop 
SS Ds eee PY = EPY (7) 
n 
n 
Therefore, the optical potential has poles at energies equal to the 
eigenvalues of Eq. (6). In the case of well-separated resonances 
one can break up the optical potential into two parts when E is near 
one of the €, say €.. Thus, Eq. (7) may be rewritten as 
0 0,|H g Gs|H 
r Hpa | n) « n| QP _pl)py- - Hoa | s) « s | QP py 
Hop K- € E-€ 
n S 
n#s (8a) 


See ey definition of H’, 


$,)(85|H 
(H’-E)PY = (ee) PY . (8) 


F-€ 
S 


The operator H’ is that part of the Hamiltonian which varies slowly 
with the energy E and gives rise to non-resonant scattering while the 
Sperapor Hpg!®s)(9s|Hyp/E-€ . pives. the rapidly changing. oriresonanw 
COucr i 2uGion tO the scattering. 

The rirst step in solving Eq. (8b) is to solve the homogeneous 


portion (H’-E)P¥ = 0, where the + indicates that the solution 
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satisfies outgoing boundary conditions and minus indicates that the 
Pelucion Satisfies incoming boundary conditions. The functions py 
describe non-resonant scattering Since they are associated with H’. 
It is assumed that Py exist and are known. Given the set of o from 
fae (6) and py one can investigate the T matrix (see Goldberger and 
Watson, 1964) and thus the scattering cross section. From the out- 
going homogeneous solution BAe a Green's function can be constructed 
and an integral equation for the total outgoing solution P¥* can be 


@erived, namely 


1 Hyg l#s)(8s1Hypl PY") 


Z 


I 8 fe B-€ 
S 


(9) 


pyt = pyt + 


Operating on the left with (Bs [Hop on both sides of Eq. (9) one finds 
the relation between (Bs | Hyp] PY") and (Bs|Hop|PYS) which, after 


substitution back into Eq. (9), allows one to write P¥t as 


Gavi gel, Hee Pye 
pyt = pyt + —1L _ ha Bs) (Bs| oP 3) (10) 
an B-(e + a’s))+ i i) 


where, using Feshbach's notation, 


S ae (is 1 
Ae) ao Ve OslHop Hpg| 9s) (11) 


Since the inverse operator 1/E+- H’ can be written as 


ae ‘seen 


ae . 2(—). imo be He ) (12) 
H 


tie explicit delanitions of Thevenerry shalt A) and the level width 


pis) 


Ale 


AG “ (Bs| Hyp P(=—s7) Hpg | 9s) es) 
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anc. 


(s) _ , / 
le = On (Bs |Hyp6(E-H ) Hpg | 4s) (14a) 
cr, equivalently, using the expansion Hpg| 9s) - fe(e/)| Pe (e’)) aK’ 


S 2 
r*) = on | (8s| Hyp PYS(Es)) | (14) 
The T matrix for the effective Schrodinrer equation (elo) wats 
(Xo{T|Xo) = (Xo{V| P¥*) (15) 


where X, lee tie Ee VeeyIUMNeULONe Ol bhai gals Of H’, Galled Hos which 
describes unperturbed motion of the incident particle and target 


fyovem, that is Ho? and where according to Eq. (8) V is given by 


; tpg! #s) (8s Hop 


(16) 


Vv 

non-resonant R-e¢ 
potential S resonant 
potential 


Inserting PYt from Eq. (10) into Eq. (15) and recalling that ese iG 
the outgoing non-resonant scattering solution of (H’- ee = 0 and 


mm@ererore Satisfies 


+ ile 
Ie are re HY Peake oa] (17) 


it is seen after some manipulation (see Goldberger and Watson, 1964, 


Sec. 5a) 


(Xo SG) = Xo = eee) 


non-resonant 


(P¥5| Hp | 8s) (Os | Hyp] PYS) 


: -(c, +a‘®)) + = ifs) — 


i 


resonant 





eal 


mince the cross section is proportional to the square of the T matrix, 
it will show a Breit-Wigner form and will change rapidly at pac + ais) 
with a helf-width of me) When E is not near €. oF ais) the non- 
resonant term will dominate and the cross section will vary smoothly 
with the energy. 

Thus with a rigorous formalism applicable to both atomic and 
molecular cases it has been shown that the eigenfunctions of the 
effective Hamiltonian H.., Eq. (6), do indeed lead to rapid changes 
of the Breit-Wigner type in the scattering cross section when the 
energy E is near an eigenvalue of Faq: It is assumed that Hog does 
have bound-state solutions in the vicinity of the energy FE. If this 
1 not the case, the resonant term will obviously vanish at that 
energy. 

In order to interpret what iS happening physically when a 
resonance occurs it 1S necessary to form a wave packet centered about 
the energy E. To determine with good precision the time at which the 
particle enters the “internal region," where Q¥ is non-zero, a sharp 
wave packet involving a broad range of energies AE must be con- 
structed with the assumption that AR >>I’. Choosing t = 0 as the time 
atv which the incident particle enters the internal region, the part of 


the packet due to the Q¥Y* part of ¥* at t = O may be written as 


av(r,t) = | c(e) avt(e) e ™* az (19) 


AE 


where C(E) is a slowly varying function of energy. From the 





Be 
expression 


QYy = ——--——— Hs pyt (20) 


(s) 
2 aay amas Ch WSC) RD (21) 
E ( +h ) += 40r ° 


To see this one has to introduce the operator Q (in Eq. (6b)) before 
the Hop operator in Eq. (20) and the expression (10) for P¥*, and 

select all resonant terms of the resulting expression which have the 
factor 1/E” - e.. In the final step one has to replace (Bs |Hyp| P¥o) 
By its absolute value (changing in this way only the phase factor of 


B.) and to realize that according to Eq. (14b) (os lipl¥s)| = Yr’) /on, 


Wierefore att =0, EWE 


L fe pls) 
2 n(s) 
Q¥(r,0) ~~ C(E) 7 ae) ae ee g (r) dE. (22) 
AE (cra?) - 2 - 5 ar 


At time t this internal wave function becomes (since v'(E) are elgen- 


al a ie sya 
(s) : rts) 
Q¥(r,t) = o(e, TA Ja.) 5) 7 5) gE 25) 
(e+ at | - fF; - 5 aa 
AK 


the assumption being made that AE is broad compared to the resonance 


fume ions of H) 


5 ae , 
at oe +a ) out small compared with the distance between successive 
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resonances. With the usual approximations for sharply peaked inte- 


grands, the integration limits may be extended to -o© and +o yielding 


QY (r,t) n= 5 ar’) o(., +a(®)) B (r) 


xr 
(2S 


we a7 tat 
4 da 
( + ae Sai 2 - Hw 
= 00 Ss Oo 


This familiar integral (in time-dependent perturbation theory) can 


be evaluated to give (see Goldberger and Watson, chapter 8) 


2 eaCs) (s) 
FG) xn if o(.+ Jo.) 
t>o 
: wi/alesta’®))t 2) jax t 
Qv(r,t) = (25) 


O oa <0) 
The Q¥ part of the wave packet thus comes into effect as the particle 
Comes into the region of the target. In this region the QY part of 
the packet and thus the Ntl electron system resembles a stationary 
eieentunetion Of H.. except for the fact that this statiomery, stare 


QQ 


is modulated and decreased in time by an exponential decay factor with 


mean lifetime T = n/p hs), It has therefore been shown that the Ntl 
particle system does exist as a quasi-stationary localized state for 
an average time t. If one looks at Messiah's development of radiation 


theory, (1963) where projection operators are used to separate the 


light field from the matter states, one sees that the use of the word 
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"stationary" in describing the states of negative ions is either 
rigorous or approximate in the same sense that one calls "stationary" 
any excited state of an atom or molecule that has a finite spontaneous 
emission probability. 

The beauty of the projection operator formalism is its 
generality. It does not depend on the choice of H, and Q and P are 
Only required to satisry P¥Y —¥Y at large distances, so that the 
formalism can be applied to any decaying state problem. Note, how- 
ever, that both the formal theory and the quasi-stationary methods 
require that the physical situation be inserted, either by chemical 
intuition or by proper choice of the projection operators P and Q. 
Mois 1s the price which must be aaa for not being able to integrate 


the Schrodinger equation exactly. 





25 


CHAPTER III 
RESONANT STATES OF HS 
iIntrodue t 1on 


In this paper the quasivariation method developed simultaneously 
by Taylor and Williams (1965, hereinafter denoted by TW) and O'Malley 
and Geltman (1965), and the stabilization method discussed in the work 
of Taylor, Nazaroff, and Golebiewski (1966, hereinafter denoted by 
TNG) are applied to the ab initio calculation of quasi-stationary 
states of Hy Results are compared quantitatively and qualitatively 
to dissociative attachment, vibrational excitation and electron 
scattering experiments. 

For the definitions of, and distinctions between the three types 
of resonance calculated here (Single Particle and Core-Excited Types I 
and II, hereinafter denoted by SP, CEl, and CE2, respectively); the 
usefulness of these definitions and a discussion and derivation of 
the variational principles used in the present work, see the paper 


of TNG. 





¥* 
This chapter is to be submitted verbatim to the Journal of 
Chemical Physics for publication under the joint authorship of 
James K. Williams and H. S. Taylor. 
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Experiments and Results 


The experiments simplest to explain with the calculations pre- 
sented here are those in which resonances are seen in the total and 
inelastic electron-hydrogen cross section. Such observations have 
been made by Kuyatt, Simpson and Mielczarek (1964, 1965b, 1966) in 
transmission experiments, by Golden and Bandel (1965) with a modified 
Ramsauer technique and by Heidemann, Kuyatt, and Chamberlain (1966) | 
in inelastic scattering experiments. These resonances are attributed 
to the formation of two short-lived states of He and the individual 
peaks in the cross section correspond to definite vibrational levels 
of the resonant molecular state. The well-developed vibrational 
structure indicates a lifetime sufficiently long for molecular 
Miprations Of the order of 10°13 eae 

As was noted in all publications on these resonances, there is a 
striking resemblance of the vibrational spacings of the two resonances 
themselves, and also of the spacings of the resonances and the spacings 


of the four H, states ol ; ay. ; cn and gy, ‘ 
2 u g u & 


sal the earliest of their publications, Kuyatt et al. report only 
One resonant state. They later refined the experiment to observe two 
close-lying states. In TW, a rough calculation was carried out which 
agreed with the one resonance observed at that time. The general 
results of TW are still valid but are superseded by the more detailed 
calculations and analysis of this work. In this work, comparisons 
will be made with the later experiments of Kuyatt et al., which differ 
somewhat in energy scale from those of Golden and Bandel. This is 
admittedly done because of better agreement of the calculations with 
the Kuyatt work. 





all 


ies INcweaves, as discussed ian ING, that the observed resonances 


are CEl, an electron in the field of a mixture of the four H, Staves 


above. Based on physical arguments, TW correctly credited the first 


2. y : : a 
resonance to a Le Seo Ia) ero) COnlasuravron oT He em gel 


electron in the field of the I state,~ 
Bites possible contreuratvions that result trom addition of an 


excited orbital to the four Hy states and could cause reSonances are 


listed in Table 1. Non-autoionizing configurations are omitted. In 
Mremcdird Section, the details of the mixing of those configurations 


with the same symmetry in both the quasivariation and stabilization 


3 


methods are reported. 


cance wewere UuSsingsopen Shell calculations, 2 prime omicedune 
midicate an orbital of slightly larger average radius than the 
unprimed orbital. The orbital notation is the standard Mulliken LCAO 
Mmeiveabion and 1s used for Simplicity of presentation. Any particular 
Mulliken orbital can be thought of as the dominant term in an 
Searanscion Of Che molecular orbitals, used in Unis paper and explained 
meme lacler section, upon a complete set of Mulliken orbitals: 


3on1y the two OplSnyepy) Mucp_}/ resonances and a /porviomec™ 
eurve D) of Fig. 1 were calculated with the quasivariation method as 
femieas the Stabilization method, Ihe latter was used for, aliaorder 
femeeeves. Ihe reason for this is the difficulty encountered am 
orthogonalization of the wavefunction of the resonance to the wave- 
functions of the states of Ho of lower energy whenever the main con- 
figuration(s) has two orbitals which are not of nx type. This 
orthogonalization is a requirement of the quasivariation method, but 
is not required when the stabilization method is used. The difficult 
orthogonality requirement will probably mean that most future cal- 
culations will be done by the stabilization method. Added to this 
difficulty in the quasivariation method is the fact that Che and SP 
resonances can only be calculated by the stabilization method, as 
puescussed inwlNG. 
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TABLE 1 


SOME POSSIBLE Hy RESONANT STATES 


Hy Parent Ht, States 

2 di 
Le eet oe! e I © qh 
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H, x2, (Kolos and Wolniewicz, 1965). 


+ 
One configuration result for H, 2 oe sone Eivas 


+ 
curve lies about 0.32 eV above the "exact" 2h curve 


(Kolos and Wolniewicz, 1965), but is used as a 


2 
PMESecoOm TOUT Aton plus can Exirasel cement. 


Comparison wich the HH. callcullauion eonstetrnceo, 


Two configuration result for H, °l (lo lx +lo ln ye 
u eee) fe 


2 

The minimum of this curve lies about 0.14 eV above the 
+ 

experimental minimum. The H, aD Curve TresmsonclLose 

CO this one thay they are indistingeulsiaete on cias 


scale. 
+ 
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The outstanding result of the present calculation for these two 
resonances was that the lower resonance is almost exclusively 
3 + oO a / @ ry a 
I e, mi ee ge geet while the next resonance is I + e, 


a yeP ea: In Table 2, the vibrational levels of these two 


eyo isn 2p 
& & u 
states are compared to the observed levels of Kuyatt et al., and good 
agreement is obtained. The calculation also demonstrated a number of 
physical effects which were discussed in TNG. First, it was found that 


even though the two CEL eI Ie states were mixed with themselves, 


+ 
with all the states of the same symmetry 2, in Table 1, and with many 


configurations of non-stable HH, 


states of lower energy, the ratios of 
the off-diagonal elements of the Hamiltonian matrix to the diagonal 


elements were O(107°) or less between the two spin doublets themselves, 


and between each of the spin doublets and the large number of other i, 


configurations. This indicates the expected extremely weak inter- 
actions (low correlation) between electrons in the compound state. 
The electron affinity of the lower CEl resonance is 0.8488 eV relative 
to its parent “I state, and the higher CEl resonance has 0.8435 eV 
binding relative to its Ml, parent, or 0.276 eV relative to the A 
CFl resonances were not attainable from configurations repre- 
eenuting an electron in the fieid of The other Hi, states in Table l. 
No extensive search was made in this work for CE2 resonances belonging 
to the states of Table 1, although to check for their existence, one 
good stabilization calculation was done at the equilibrium internuclear 
separation of the four parent states (R = 1.95 au, 1 au = .52917 A), 


for the configuration an electron in the field of on “Dto,1s0,260,28/. 
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TABLE 2 


COMPARISON OF EXPERIMENTAL 
AND CALCULATED VIBRATIONAL LEVELS 
(energies in eV) 


Lower Resonance Upper Resonance 


Cake, Obs. Cale 
hewele “sp leueane bevel oplitting Slevel ssp iser rad 





als GRO eae te Ibe Ss 
28 . 20 .26 . 29 

ele 56 eee ite fe IS 
28 29 o2g 28 

11. 84 11.66 11.99 12208 
ee ey 28 28 

ee ily oe. ave7 12 eh 
. 26 .26 . 26 mer 

en 37 12.19 1252 Ne. Sts 
5 peu an 26 

12.62 12.43 Heal 12.84 
24 722 520 (25 


12.86 12563 Leo 13.09 
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This root stabilized at an energy of -.72911 au (1 au = 27.21 ev), 
a 
0.24 eV above 3 -3 eV below Is . 13 eV above De Aid (lees sey: 


above the ground vibrational level of H roy, the reference zero level 


2 


fmreall the calculated and experimental results. This root strongly 
madicates the existence of a CkK2 resonance with this configuration 
which parallels its on parent with a slightly higher energy in the 
Franck-Condon region. 

Again this eee configuration does not mix with any of 
the ee ee ee types and the resonance is expected to be quite 
broad relative to the CEl resonances as discussed in TNG. It is likely 
that this resonance will only be observable in inelastic scattering 
experiments leading to its parent, however it is possible that thres- 
hold structure will obscure the observation. Because the two CH1l 


resonances are made from excited H,, orbitals, it is to be expected 


that these resonant states are more strongly coupled to excited state 


exit channels of Hs and this accounts for the enhanced signal due to 


these resonances in elastic scattering (see Heidemann, Kuyatt, and 


m 
Chamberlain, 1966). Z 


ly 


In the next section, an even more tentative CE? resonance, with 

the configuration of an electron in the field of EL, “di dgls 0g2sdy2s 
S 

is reported. The uncertainty of this resonance causes us to discuss 
Mimonly Sin tne Next SeCcGlon, bul 1bsis ment 1lonedehere inet he scat 
tion that many readers may not be interested in the detailed dis- 
cussion of the calculations. Work is now in progress on this and 
Other CE2 resonances expected to lie in this energy region. We 
emphasize that the two CEl resonances and the single point at 11.85 ev 
are definitively stabilized and therefore represent final results. 





3h, 


In Fig. 2 the two CEl resonances and the four aforementioned 
states of H, are shown, as well as the single point for the CE2 


resonance (denoted by a cross). The solid lines represent the H CEl 


resonances, the dashed lines the i, states, and the dotted line 
through the cross is the postulated potential curve of the Ck 
resonance. 

ie next Gype Cf Experiment to be discusséd is dissociative 
attachment. The measured cross sections for the formation of the H™ 


ion in the process 


H +H 


e +H, > (HZ) > 
eS 2 H* 4 H- 


are given in Figs. 3 and 4. In an as yet unpublished private com- 
munication, Dowell and Sharp have reported finding additional structure 
in the 8-12 eV peak. Besides the general Gaussian shape and the dip 

in the cross section at 11.2 eV, they have found a new series of dips 
on the high energy side of the peak, which have the feature of being 
spaced analogously to the spacings of the vibrational levels Kuyatt 

et al. found in their lower CEl resonance. In TNG, it was shown how 
these experiments could be used to completely and definitively predict 
the position of the resonance potential curve causing such structural 
features as contained in the 8-12 eV peak. First, the general Gaussian 
shape, kinetic energy analysis, and isotope effect could be caused, as 
a result of Franck-Condon factors, by any repulsive resonant state 

With a potential curve which enters the Franck-Condon region of the 
ground vibrational level of the H, ground state at about 8 eV on the 


2 
right, leaves it at about 13 eV on the left, and goes to the known 





Fig. 2.--Comparison of the upper core-excited resonances to 
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Fig. 3.--The dissociative attachment cross section of 
Schulz (1964). 
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Fig. 4.--The dissociative attachment cross section of Rapp, Sharp, 
and Briglia (1965). 
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stable separated atoms H and H™ at an energy of 3.75 eV. The dip at 
11.2 eV was most easily explained by assuming that this repulsive 
resonance potential curve could cross the repulsive H,, De state at 
ae. 2 eVv,? thereby changing from a CEKl to a CHe@ resonance, and opening 
a new exit channel H + H + e to the system. Since, from the work of 
Kolos and Wolniewicz (1965), it is known that the H, i state has a 
relative energy of 11.2 eV at R~ 1.29 au, it can be concluded that the 
Hy, has the same energy at this ae The opening of the new exit channel 
causes more of the H, resonances formed at this energy to decay in 
channels other than that measured by dissociative attachment, and there- 
fore causes a decrease in the cross section exactly at 11.2 eV. The 
reason that the cross section recovers its aehesaet envelope is illus- 
trated in Fig. 5, where it is seen that only at the crossing point are 
the Franck-Condon factors of sufficient magnitude to make possible 
eurve jumping to the os ie ee Mey ia rie lies 

The Dowell and Sharp result could most simply be explained by 
assuming that the resonance curve proceeds upward in energy until it 
approaches the lower “Dro, 1sx, 2pm, 2p/ CEl resonance of Kuyatt et al. 
tangentially from the left. As shown in Fig. 6, if this were to occur, 


i: : : vias 
Since both states are “a one could get an inverse pre-dissociation 


from the repulsive A state to the bound A, state at energies equal to 


ee nedis in order to have the maxima of the nuclear wave 
functions coincide, the crossing would occur at a slightly lower 
energy, as shown in Fig. 5. As a consequence, the R value for the 
crossing point predicted by the dip will be slightly smaller than the 
actual R value at crossing. 
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Fig. 5.--Comparison of Franck-Condon Factors. 
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Fig. 6.--Schematic diagram for inverse 
pre-dissociation between resonances. 
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the vibrational levels of the upper resonance. Because of Franck- 
Condon factors (see Fig. 6) the postulated inverse pre-dissociation : 
could only occur at vibrational levels of the upper resonance, and this 
Sema cause exactly the dips in the H’ formation cross section observed | 
by Dowell and Sharp. 

since from experiment one knows the vibrational constants of the 
Kuyatt state, its potential curve can be determined empirically and 
miesviprational levels drawn in. By noting the energies and positions 
weewich the vibrational levels and the empirical potential curve 
iawersect, and by ignoring those points which do not correspond to 
f@ipseein the Dowell and Sharp experiment, one is Petuewa via see on 
points (denoted by crosses in Fig. 6) through which the repulsive 
resonance potential curve must pass. If one now smoothly connects the 


"cross" points of Fig. 6; the point on the H mut curve; the point at 


a 
which the 8 eV constant energy line intersects the large R, vertical 
Franck-Condon line; and the point at infinite R for the known H + H™ 
energy, a curve (such as curve D of Fig. 1) is obtained with little 
ambiguity. In TNG, a conglomeration of tentative calculational, group 
theoretical, curve crossing and physical arguments were given to pre- 
dict the curve, which, it was then pointed out, could explain the 
experiments. The opposite point of view is now being stressed, that 
once the type of phenomena that can occur in quasi-stationary molecular 
resonances are understood, and these are quite closely analogous to 


normal spectroscopic phenomena (i.e., dissociation, inverse pre-dis- 


sociation, pre-dissociation, etc.), one can predict where the resonant 
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Beate lies completely from experiment. 

As also explained in TNG, simple Mulliken type arguments can be 
memoaned with the core-excited model to predict a configuration for 
ime resonance of “Yio 1so,1s0, 1s’. In this work, using the stabili- 
zation method, the energy of such a configuration was calculated, and 
feeeeresults which are presented in Fig. 1, curve D and in Table 3 are 
[memecrirec. agreement with the curve one would anticipate from the 
application of spectroscopic ideas to the experimental results. As 
will be shown in detail in the following section, this state also had 
maeremely low correlation, since it did not mix appreciably with any 
configurations which would normally be regarded as correlating con- 
meourations in the configuration interaction picture. Some mixing did 
@eeir with the stabilizing configurations, particularly in the CE2 
pemcton of the curve, imdacating that the CEe will have a larger 
width (be less stable) than the CEl part. Again it is stressed that 
titeis the lack of correlation in the resonance complex that enables 
one wavefunctions determined by quasi-stationary methods to give such 
Sood agreement for the energies at which resonant phenomena occur. 

Work as presently in progress which willeuse the waverunericns 
Seucrmined in this work for the “Dto Iso, 180, 1s/ state in conjunction 
with the formal theory of dissociative attachment and vibrational 
excitation as developed by O'Malley (1966) and Bardsley, Herzenberg, 
and Mandl (1966a,b) in the hope that this will lead to quantitative 
agreement with experiment for dissociative attachment and will explain 


why  oChULZ US Unables to Observe any evidence or the “Eto 1s0,1s0,18/ 
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resonance in electron impact experiments where the cross section for 
excltacion from the ground vibrational level of H,, HS) (en.dealicie(ol 
mepravional levels is measured. Since the resonance potential curve 
gees Cross and go above the =e channel in the Franck-Condon region, 
ena tous lies very close to the H, state in the center of the Franck- 
Senden region, it is unlikely in any case that transitions throtigh this 
mesonance would make significant contributions to the vibrational 
=aetuablon Cross Section to the lower i, Vibrational levels. 

ine shawe = teoLepe Cliecr and kwnebic enersy -analcieson ceme 
14.2 eV peak were explained in detail in TNG, and credited to 
means’ t1on to the dissociation limit of the two CEI “Teo 1sm, 2pn 2p! 
resonances discussed above. The calculations presented here reinforce 
emerare in complete agreement with the ING arpumentcs hence the argu- 
fees will not be repeated here. In .a future paper, the wavefunction. 
wakewlacred [Or Liése@stales as a funculon Of internuclear distance... 
also be used in conjunction with the formal theories of dissociative 
attachment mentioned above. 

It was postulated in TNG that the vibrational excitation results 
of Schulz (1964) and the dissociative attachment cross section, kinetic 
energy distribution, and sudden onset for the peak occurring at 
3.75 eV (Schulz and Asundi, 1966) could be explained by the existence 


of what must be a quite broad SP resonance with the configuration 


2 
2, 6 lso 1s/o Ps, 1a Ofle electron an tue freld oF x xiyt 6 
ug g u u fis g 





ene reader is reterred to ENG for asdsscussion Of own this secon— 
figuration was chosen; why it is expected to be broad; why the reson- 
ance is difficult to observe in transmission, mocified Ramsauer, and 
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The stabilization method has been applied for this state. A 
good ground state i, FUNCT LON was ~enOsen,) 176. ,.c three cout ieurarion, 
well correlated function which gave an energy at the minimum of the 
potential curve of -31.76 eV as compared to the experimental energy 01, 
~-31.975 eV. The three configurations used were Some 3 oe liscier 
and n 2p," 2p’, - The first configuration gives the main contribution 
around equilibrium nuclear separation Ro: The second is an axially 
Serrelating configuration near Ro and ensures proper breakup into two 
Meacoms as Ro, where 1t becomes as important as the Neral con- 
meouration. The npn 2p/ 1s an angular correlation eonmieur at ren. 
Sitemsuperscript c is added to these orbitals to indicate that they 
meve une same average radial extent as the ora orbital, and are 
eeecilally Optimized to produce the best i, SYOUNG Sale, as OppeCseemme 
the gis and 1 ep orbitals of the previous discussion which are 
Similar to large excited orbitals of H,. 


a 

Ror the ae i, Calculations ape ire -cucernuen ois was added to 
meeecrconripurations to give three three—-electron coniigurations of 
symetry “2”. This was then stabilized as discussed in the next 
section and resulted in the curve E of Fig. 1, as shown in detail in 
Bee. * and tabulated in Jable 3. it 1S notable tray even though one 
needs at least three configurations to represent this resonance well 
im 15 Only cecause one must represent the target H, state to a high 


@e ree ot accuracy. This type of configuration mixing does ner 


Cetal Cross Secu@on Experiments, and whne=seneral physics of the 
Situation. 
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Fig. (.--The single particle H 
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Mmiaveate Strong correlation between the new electron and the original 
meeelectrons. On the other hand, it was found that some configuration 
teeing was necessary before the root stabilized. This situation is 
peaeerous to that of the CEe portion of the “Yo, 1s0,,180,,15/ 
mepulsive resonance, and again the extra mixing with the stabilizing 
Pemreurations indicates thal this resonance is guite broad. In any 
case, one finds a wavefunction with a small number of configurations 
whose eigenvalue was unaffected by further variations of the wave- 
mum@erion. How much of this configuration interaction is due to 
fee, and how much 1s due to a limited orbital functional form is 
ampossible to tell. 
Moe iiteereceime Teabures of this resonance are; 

i) That at large R it connects smoothly to the H + H™ 
a uaee Ccaleculaced by Taylor and Harris (1963) uSing a true vari- 
ation method since at large RH +H is the bound, lowest state of 
on symmetry. 

ii) That the potential curve for this state, once inside 


al 
the curve for H, X ae essentially parallel sence wou Csi) Cilmi eve. 


e a 


and turns up sharply at R~1.0 au. Since its shape is like that of 
He, its harmonic vibrational constant will be similar to Hy» and thus 


the lower vibrational levels will be at ~ 1.5 and 2.0 eV above the 


ground vibrational level of H, X ae 
The calculated state could well account for the dissociative 


attachment result of Schulz and Asundi (1966), but it should be 


pointed out that the dissociative attachment experiment in this case 
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feres nO information about the position of the potential curve for the 


mesonance inside the target H. potential curve. 


2 


From the vibrational excitation experiment, Schulz (1964) 
Peevulated the existence of a very broad resonance about 2.3 eV above 


1 
mae ground level of Hy X ue The VSesiite presented Here ts ecer rcmml: | 


lower than this. The only significant test of the calculated state 
Haat cOuld be carried out, but for physical reasons, one which might 
not be extremely sensitive to the detailed form of the wavefunction, 
is to put the wavefunction calculated here into the formulae for the 
feoravional excitation cross section as presented by Bardsley, 
Herzenberg, and Mandl (1966a,b). This has been done by the latter 


authors using a wavefunction of significantly poorer quality than the 


vf 


Mmaemorescitea here and nevertheless resulted in a calculated Cross 


feellon which 1S in quite acceptable agreement with the observed one. 


‘It demonstrate this point, one need only note that the 
Herzenberg et al. wavefunction gives rise to an Hs energy of -30.44 eV 
and an equilibrium internuclear separation of 1.68 au. The wavefunc- 
tion used in this work, gives an energy of -31.76 eV and an equili- 
brium separation of 1.4 au. The “experimental” values are -31.975 eV 
at 1.4011 au (Kolos and Wolniewicz, 1965). Because of the different 
qualities of the wavefunction and the wide variations in Rg, it is 
meaningless to try to compare the results of the present HS SP calcu- 
leblon With that of Herzerberg et al. Although one does net have a 
Warlavilon principle which says the lower the energy the better the 
iinccvon in these calculations, 1b 15 reasonable to assume that the 
more extensive and better functions used here will lead to more reliable 
meoults. The agreement with experiment in the Kuyatt et al: resonances, 
and the ability to explain the fine structure due to the repulsive 
Pocomeance. which as completely Macking in “tie Herzencera et “al. results, 
Sec Lo substantiate these claims. 
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Bie oiis laboratory the Herzenberg et al. calculations for the vibra- 
tional excitation cross section are being repeated with improved wave 
me ctons to see how sensitive the calculation is to the quality of 
the function. 

The qualitative comparison of this state to the energy of Schulz 
Semeeave difficult since the SP state is so broad. in fact, in the 
feoralional excitation experiments, since the resonance is so broad, 
fieenuclei do not have time to make definitive molecular vibrations, 
and Schulz was not able to distinguish any meaningful energies to 
determine the empirical potential curve for this stake. 

When the a CH, X 2) ois root was being stabilized, it was 
noted that a second root stabilized equally as well. The curve for 
mars is fF’ in Fig. 1, and the upper Hy curve of Fig. 7. It seems to 
dissociate to H + H’* where H™™ is probably the very broad H™ resonance 
1s 2s calculated by Herzenberg and Mandl (1963) to have an energy 
between 2.0 and 4.0 eV and a width of about 2 eV. This i Stale, vacam., 
moll be very broad, and cannot be seen in a dissociative attachment 
Paperiment since a long-lived negative ion is necessary. As discussed 
Meee lecause of its breadth at also might be expremely difficult to 
@etect in the total scattering cross section. Since it oceurs at about 
[ee mms could possibly contribute to Vibrational excitation.  Onesmient 
speculate that both the 1.2 eV and 3.0 eV SP resonances contriosute to 
give the broad vibrational excitation peak estimated by Schulz to be 
meeeered ab 2.3 €V.° Nothing détinivive can be until calculations or 


widths and the vibrational excitation cross section have »een made. 
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Pa IMceresuine scr Lect Observed in these SP calcubations was that 
mierwextra electron of the higher stable root had a smaller effective 
orbital radius than that of the lower root at small R Co euls Durer 
fee KR the situation became reversed. 

When stabilizing the repulsive Yo 180 180 1s/ Stace, id aime 
MmenlLavions reported here, because of the open shell methods in which 
imeeesoin doublets are used, a configuration 2) 6 1s0,,1so,,18/ (electron 
mamoie field of BD) was automatically included. During the stabili- 
mieion process, a second root corresponding to this configuration also 
meebilized, and it was clearly a broad, CE2 resonance. This resonance, 
fem tne other CHl and CE2 cases, was well stabilized and highly 
Pereorrelated and fits the rules given in TW and TNG for a resonance 
femete t1eld of the B state. It should be difficult to observe for a 
Mieecr OL reasons. As discussed in TING, the breadth of this resonance 
fovers its being looked for in inelastic experiments to vibrational 
MwevelsS Of its parent B state. However, besides its breadth, the Farce 
Re of ~ 2.45 au and the small w, of ~ 0.13 eV indicate that the 
mesonance potential curve will be in such a@ position relative vo the 
Becential curve of the H, ground state that Franek=Condon facters will 
allow transitions only to the highly excited vibrational levels of the 
resonance (v ~14) and any structure may be obscured by the closeness 
of the spacings. It is impossible to tell from the published curves 
of Heideman, et al. (1966) whether there is any structure just above 
Eimer excitawron tmreshold, and Since the wide ais neu kom, it is 


impossible to really argue whether or not such structure should indeed 


be observable. 
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Wiewctrabilizaplon Chlunissroou, which fits the Ch model ef TW 
and TNG, brings up a very interesting and unanswered question; it 
mines tnat in the Gase of@eore-exciled resonances a well stabilized 
root appears fastest whenever a configuration is used which puts the 
extra electron into the same type of orbital as the excited i, orbital 
fe.e., 0,(0,,)" and 3, (m,)°]. Othermsconiirurations do not svapillize 
rapidly. Why this is so, we do not know, and it is freely admitted 
[mreemetne extent of calculational experience so far is not sufficient 
Pemsay definitely that these other configurations do not stabilize, 
mi there are strong indications that this is the case. Another 
merange effect arises when one adds many configurations solely for 
the purpose of testing the stability of a particular roou. lf the 
root under investigation is one which is truly stable, two or three 
meapilizing configurations are usually sufficient. On the other hand, 
@paer roots of the secular equation do eventually become constant 
after large numbers of configurations have been added, simply because 
iMo@eeeconiiguratlions made of orbitals which can mix with the configura-— 
tions corresponding to these roots have been exhausted. A subject of 
imine research, especially when reliable widths cangbe produced, wil 
Sem cuudy Of the significance of these more slowly ssca0ilizingesvauvee. 
ijesevunems stress that the states reported ii tie presentawor, 
womverzea rapidly, and tie Const onme.1SlS "onl eobou. Griereless 
Beaple roots which were noted in thescalculations. 

In the tables and graphs the energies are reported in electron 
volms relative to the experimental ground vibravional level of i, in 


the ground state. When enerries of resonances relative to energies of 
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H, eee reported, the H, energies are those obtained with the same 
wave function as for i, with the exception that the linear coeffi- 
ments are re-Optimized and the third electron is allowed to go to 
Mminity. This is the only way to get a true comparison and does not 
yield misleading results since the wavefunctions used in the i, case 
eave well over 99% of the total energy and an appreciable fraction of 
ie correlation energy. In view of the experimental errors reported, 
imeroving the calculation of resonance and target energies to a degree 
Saesophnistication comparable to that of Kolos and Wolniewicz GICG Sn 
fem do little to improve our understanding of the basic physics, 

and therefore would not justify the many man-hours and computer 
emecises involved. 

A point of great interest arises when discussing the non-crossing 
rule of quantum chemistry. As noted by Bardsley et al. (1%6), and 
“ne present authors, many unstable potential curves cross the stable 
mevential curves. This would seem inadmissable by group theoretical 
Peeuments which say that the correct potential curves must be drawn by 
Pennecting the aa lowest root at all R. These arguments have meaning 
in the variation method, but not, however, in the quasi-stationary 
methods used here. Bardsley et al. show this by reference to the 
work of Mandl (1966) in which it was proven that for resonant states 
calculated by the quasi-stationary complex energy Siegert method, the 
real part of the energy can cross (and this is what is calculated heme) 
but that there is no crossing in the complex plane. While the authors 
have no argument with this correct interpretation, we prefer to pre- 


sent a more physically based model. To be specific, consider the 
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2b 9,180, 1s018/ repulsive reésolany stave, which avy larce R dissociates 
womaeand H'. This separated atom pair in turn lie about 1.0 eV below 
a pair H and H~*, where the H-* represents a SP resonance of H- dis- 


eiuesed by Herzenberg et al. The upper pair also has the symmetry 


+ 
Sy Meo 18 06 1so 1s’o 2s. 
& & & E 


In the united atom limit, the “Dto 1so1s0,18/ state goes to an 
obviously core excited state of He”, while the H and H~* pair goes 
to Heo “Sle es. Ihe latver 1s clearly a very broad SP stave o1 Her 
and lies significantly lower than the core-excited state (about 19 eV). 
This means that the rapidly stabilizing eee ie ee root and the 
unstable eS rOOs, whieh are mixed GO teste the sravilizaciom 
mapcne former, should cross av a reasonably large internuclear distance. 
itars is exactly what was observed in the calculation. The solution to 
mimes solLlLémma is to realize that curve crossing is a dynamical as well 
fea rroup theoretical problem, and that the nucléar motion will thus 
follow the diabatic curve (wave function envelope) as opposed to the 
adiabatic curve (energy envelope) whenever the velocity of the nuclei 
is significant and the electronic interactions between the two wave- 
mmerlons corresponding to the crossing curves is insieniticans. In 
fm@escase Of Core-excited resonances, it has already been polnved out 
peveral Limes that the mixing between the stable and added stabilizing 
SOntigurations is 5-10 orders of magnitude smaller than the diagonal 
energy elements. (As discussed in TNG this fact is physically obvious 
and fundamentally the reason for the success of the quasi-stationary 


methods.) This mears that in the problems considered here, the 
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mr peuic curves are the ones which must be compared to experiment. ‘The 
semeral procedure, themetore, tseto follow the dtabatlic cubwesswhich 
Can easily be chosen by studying the smooth change of wavefunction as 
Memertton Of R, and then by looking to see if there 1s any Sip- 
Mificant mixing when two diabatic curves cross. A note of caution is 
mimorder —— in processes like dissociative recombination and 
associative ionization, which are basically resonance processes, but 
where the relative kinetic energy of the nuclei is never more than 
~.03 eV at 300°K, one should follow the adiabatic curves. (A 
mmielitication 18 introduced in these problems, since the SP states 
are no longer resonances, but singly excited states of the system 
formed when an electron 1s addéd to the sooee oHisie bl Iestie akovel = yal! 
maerefore amenable to the quasivariation or variation methods.) In 
Bomaerline cases, one will have to either refer to the theory of curve 
Grossing, use Mandi's complex energy theorems, or hopefully, resolve 
the question unambiguously by comparison to experiment (one of the 
Sioces, diabatic or adiabatic, should give consnicuously poor amrec- 
ment ) . 

In this section any discussion of widths has been knowingly 
Omtced. For the aN Ved..16/0 1s) oPeand ae lso lso ls’ repulsive 

be h U gg u U 

Suaee, the widths have been@ealeulaved=by Herzenversa et al., sour ior 


meacous Ghat will pe distussed in a fugure puslicanmon, thesageunors 
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Beserve judgment as to their validity. 

itteas Clear foreG@hiresenances that narrowness of widGhe greatly 
mevors the simplification introduced by usirg quasi-stationary methods, 
Memerp sed to pure *icse-coupled scattering calculations, when one only 
fereres to have infcrmation about the resonances’ energy, width, and 
assignments. For CE2 and SP, which are quite broad {lifetimes of 
6 (107!) SEC. or less {7 Some authors would prefer to trea tie ihe= 
gum aASSOciatéd with these states “y pure scattering methods, 
eeeocially the method of distorted waves. Herzenberg et al., i:. their 
Bmeeulation of the vibrational excitation cross sectior, stress that 
they obtain results with quasi-stationary methods for SP RONG SS — 
Seeneare in as food apreement with experiment as the result of 


Takayanasi (1965) which was achieved with a distorted wave method. 


Sohne widths of resonances as calculated by complex energy methods 
~melar tO those used oy Herzenberg et al. will oe discussed in a 
puolacaLion in préparation by H. Taylor and A. Golebiewski. It will 
memenown there that the calculation of widths gpives ambiguous results 
when the boundary enclosure (see Herzenberg et als) is varied unless 
extremely large basis sets are used in the inside region. No calcu- 
iiiens have ever peen done using the complex energy method for widths 
Memvigvel extensive basis sets were used. Such caleulavions re in 
Peeetess at this laboratory, and will be the subject of future 
communications. 


Jeareie Particle resonances are expected to be even broader 
Poeauce Ol Lhe relatively lower polarizability of the ground state. 


Venous authors use the following names for what are denoted 
here by CE2 and SP: Shape resonances, virtual resonances, potential 
Beseances, direct resonances, etc. 
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pee Demuth oul wy werZenwere EL al, this indicares thacge sie two 
meeernave models “approach the true situation, that of a distorted 
m@recy and an interacting passing electron, from two different 
mimating points of view. 

The resonance formulation envisions a quasi-stationary state 
mitch must be coupled with outgoing waves, while the distorted wave 
method envisions a polarized target and a distorted ornit for the 
moscing electron. Both of these approximations, if refined, will of 
meurse five the same result. 

One can certainly say that the resonance formulation allows one 
ifemearry over the chemical concepts that are well known in spectro- 
ScOpy and quantum chemistry. The Meer freely admit that in mole- 
euler resonances which live Wouee sec or lless,eunere 15) pernaps a 
@ianicitatlive danger in drawing potential curves which depend on Strict 
adherence to the Born-Oppenheimer approximation. For such short life- 
fmmes it is dangerous to talk about potential curves, vibrational 
mewcls, 6€tc. On the other hand, as a first approximavion i guanri- 
tative calculations (as shown by Herzenberg et al.), and in quali- 
tative explanations and predictions of experiment as shown avoove and 
imecNG, the usefullnéss of the quasi-stationary approach and chemical 


intuitive ideas cannot be denied. 
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Calculations 


Calculations were done on CEKl, CE2, and SP resonances. For CEl, 
the quasivariation method of TW is adequate, as evidenced by the 
results published in TW. However, to calculate the CE2 and SP 
resonances, the stabilization method must be used. Hence a basis set 
for use in the stabilization method must be constructed which is 
physically reasonable and insures convergence (stabilization) with as 
few functions as possible. 

The calculations were done with the diatomic molecule program of 
Harris (1960) and Taylor and Harris (1963, series), modified to 
include the quasivariation and stabilization methods. In this program, 
trial functions are composed of linear combinations of products of one 
electron orbitals, so the one electron orbital basis is what is 
sought for the stabilization method, as used here. 

For a homonuclear diatomic molecule, suitable molecular orbitals 


are defined as 


Q 


»(5,asnmvp) = exp[-8&+ivo) poe een yas ee n™ {expt -on] 


i (1) 
4 (=i) “Expl © nis 


where €, n, are the prolate spheroidal coordinates; 6, Q@ are non- 


linear variable parameters; n, m are integer exponents; v is the 


ee 





De, 


meanicoal quantum number; p defines perade or ungerade symmetry of the 
meer tai>: r indexes the root of the chosen symmetry eae most cor= 
fesponds to the particular orbital B. The exponent k 1s Weed" to 
produce the g or u symmetry defined by p through the relation 
kK=-m+p +t lv | wnere p is even for g and odd fori. A complete dis- 
eussion may be found in the theses of either of the authors. 

A more compact notation 


rh. = 6,,(,a3nmvp) (2) 


Mmenow possiple, and will be used in the following discussion. As 
esample, 


ik a 8, (8,01;0000), Fea B, (,@30001), andl = 8, (5,@;00%11). 


Note that (1) is a highly flexible function with two variable 
non-linear parameters. 

One possible choice of basis®woulld’ be a Ssetroieerpi tals las. near, 
Hs orbitals as possible, and the functions (1) can be selected to give 
mepercgcllent representation of these orbitals. Por example, au hy 250 
eel au = . 52917 R), Lhe equilibrium separacvionewore tne sround eran 
cna HR, al oe orbital of the type (1) gives an energy only 0.00002 Har- 
tree above the exact energy as determined by Bates, Ledsham, and 
Stewart (1953). A basis of this type would be ideal for electron-H, 
Eeatvering calculations, but for the target H, and caleuilerions@on, CE 
Pesonances a Slight modification causes faster stabilization. The 
vasis was actually chosen by the following steps: 


in for a particular, internuclear separation &, the best 


i oe HS Orbrral 1s) 1ound As ds Known this fives an 





60 


excellent tepresentation of the innermost orbital of the 


excited H, states which iS deSirable Since core-excited 


resorances are in the field of the excited H, states are to 
pe studieu. 


Ao fe ee Ori eval te ee ene duded and. (ie Udine on fixed, 2s 


~ 


minimized to make the Dest one term fround state of H,. 


US vow clvcetaed HS GOteivais bhis One confisumagion H, result 
ISwiol the sest thay could Merachieved withsane daavenue 
BroOcran, 


Ore al. om and l sue Opbivals are 1Tound by MmanimiZaulen Linge yoe er cid 


2 
of tie 1 oe HS Orbital te make syhe Destaoncmeerm ole and 


3515 states of H,. 
u a 
4%, Higher orbitals 2 we cea ed, 30, *** 3 and 
a Ae) 3 Meme Ae CONStUrUGved Dy minimizacion Of nwo: 


rouus Of UUme seculay equaeiom vich the titel fiacmron 
ie 
orlan: Ne) yc ralo en, 
Eee je a Oe 


heldine “temic. ah) Olde dS wie coe 
The seus oO.) Orpeiecalsaonr a2 ae Symmetry fer Various values 
of R which were used in the CE calculations are listed in Tables 4-6. 
Products Of Chireevet these Orplttals are ie cdus ee Ost oer mun c 
H ealculations, and in the major cConiieurations weuresenting the Cul 
resonances the parameters of the outermost orbital were also varied 
to minimize the energies of the roots which were being stabilized. As 


all resonance energies (even those obtained by the quasivariation 


method) were checked for stabilization, this introduced no difficulties. 
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TABLE 6 


x BASIS ORBITALS* 
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Mm PP PP PPP PB 
NO 


Ee 
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SETOROREGEL. cee ERA cae Oo moea ene ol, 


b No.-nmv, e 

In -001 In /-001 en -201 3n -121 

5 Ol 5 Ol 5 (a=0 oO (O=0 
$20985.420207 » 24643, .25342 . 39992 030842 
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. 62398, . 52846 .48920, .49619 . 19724 . 086289 
66240, . 56064 ip) e)p) Siay-lele.: . 84183 08868 
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. 81159, . 68984 SO4243 7.64516 IF enuoreul . 08804 3 


Orbpltals sare madesby taking thesr -Counverpart sol 
EHese, 


u 
and the In’ to Tor the © I Stave, 


- AS minimized in quaSivariation method. 


- nmv changed to 42-1 for all R values below here. 


Sizenge Ol (Hx, 
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No parameters were varied in CE2 or SP calculations once the basis 
was obtained, as this could lead to disastrous results. 


The spin eigenfunctions used were the two three-electron doublets 


+ (2008 - afa - pac) 


‘6 


and (2) 


1 
—(apa - Baa) 
V2 


which can be seen to be C, an-BpB + C 


7 ap + Ba)-a and C, (ap - Ba)°a 


al 
respectively, or therefore one doublet composed of the two electron 
triplet plus an electron, and one composed of the singlet plus an 
electron. In the complete open shell treatment (different orbitals 
for different electrons), both spin doublets are used, and the coef- 
ficients corresponding to a particular combination of a spatial con- 
figuration (product of three one-electron molecular orbitals) and 

spin doublet gives an indication of the two electron target state To 
which the third electron became attached. 

The first problem considered was the attempt to improve the 
mesults of TW, and to see if two CE] resonances could be found, rather 
than the one reported in TW. 

Hence the first calculational step was to do a calculation in 
which pall ee nt and — ale were mixed and the ground state of H, 
projected out in accordance with the quasivariation method as dis- 


cussed in TW. The singlet and triplet were automatically mixed by 


the use of both doublets, so that only two spatial configurations 


were necessary. 
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lie resmiit wasmunal even though the two electron Spatial con- 
mecaations Of che two triplet states were quite close in energy, the 
maree electron states did not mix to six orders of magnitude, and the 
meeclet and triplet parent spin doublets also did not mix. The two- 
electron II states have considerably more electron affinity than the s 


eae 


Boacves, and both CHi resonances can thus be formed from ae plus an 
electron, with the ae plus au ebcetron yielding Cho resoilances, 
Me results, presenved resactually two-Conriguratiom Pes inc. 
3 


iL ; 
mDecause the parent target states I and a were determined with two 


Bent isurations, 4G, Met 7 ta to te In the singlet Tl state, the 


2 


4 


orbitals ln’, Ti of Table 6 were used instead of ee C, was 
~.99 and C, ~ .O2 at all internuclear distances. For the three 
eeectron calculation, the oe orbital was added, and its non-linear 
Peeame pers varied tO minimize the energy. The linear coefficients were 
mmcnanged. 

The agreement with experiment demonstrates the efficacy of the 
one-electron model (product of one-electron functions with small 


number of configurations). At several R values, additional configur- 


ations were added, including higher in Orbitals, oUt Nor additional 





Sie ouanan. in a two spatial configuration three electron 
peuple calculation there are four roous, and tol Minear coctiicienice, 
pace the roots corresponding to the triplet and Ssineler Spin parents 
do not mix, equivalent results are obtained oy using only the doublet 
1/V6 (200p-opa- pac) with the lower CEl resonance, and the doublet 
1/V2 (apa-poa) with the-upper, (Tien the linear coer ficients ane 
unchanged from the two electron calculations. This provides a useful 
simplification if the wave functions are to be used for further 
Cea cullbacvdons, Such as widths, etc. 
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binding to the TL, states was obtained. The energy did improve, but 
only to the extent that the energy of the two-electron states them- 
selves improved. In the interest of computer time, and since the 
necessary quantitative information for comparison to experiment has 
been obtained, the two configuration results are considered adequate. 
Stabilization was checked by adding configurations of the type 
eH, (x Dt) De, i.e. {C\o,0,7 4 C50,,9," zt Conn,’ } noe ny = 2, sere. 
The roots corresponding to the CEl resonances were immediately stable 
to six figures, the matrix elements between the CEl resonances and the 
stabilizing configurations were completely negligible (19078 or less) 
and the corresponding coefficients in the eigenvector matrix therefore 
were also negligible. Addition of more gs orbitals only served to 


meenuly lower the lowest root to an energy nearer the correct enersy 


Seecne pround state of H 


p>? and inserted unstable roots between the 


lowest root and the root corresponding to the CEl resonance (plus 
meaningless roots above the expected range of the CE2 resonances). 

The next step was a quasivariation calculation of the resonance 
corresponding to e3Yt for comparison to the dissociavive attacmmenr 
experiment . 

At large R values, no diffireullies were observed, bubeas RK was 
@eerecased into the Franck-Condon resion, thevamounc of \bindane for the 
third electron also decreased. Somewhere between 1.4 and 1.3 au (R) 
the binding for the third electron became negligible (0.07 eV or less), 
oad-as KR was decreased further, the potential curves for the Cll “tt 


resonances and the ie repulsive state of i, were coincident for ali 
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Heaetical purposes. This result 1s consistent with the model of 
erossing to form a CE2 resonance in the Franck-Condon region. The 
erabilization method was applied next. When the roots for K greater 
mien Or equal 1.5 au were checked, no change was observed, and the 
fmexangs of the stable resonance configurations and the stabilizing 
Seontipgurations was negligible. 

However, for smaller R values, the mixing 1s not negligible, and 
the roots corresponding mainly to the resonance configurations began 
to move to higher energy values. Stabilizing configurations of ne 
type eh ; aoe were used, for n from two to seven. “he Trocuc cormes— 
ponding to the resonance usually stabilized once the three stabilizing 


1 3 


three additional configurations were added at each R as a precautionary 


iL 
configurations Le {c a 1 aes ae io had been added, however 


check. 


ihe final result obtained was that ter tne smallest H cenciderca: 


2 
that the resonance changed from CEl to CE2 between 1.4 and 1.5 au. 


iO) au, the oa resonances was 0.2 eV above the De Sbave Of iene 


In these calculations, simple products of one-electron orbitals 

were again used, with the 3st of lo lo and ae of lo lo 60 from 
U g ou g ro) ULE UE 

Tables 4-6. For R > 1.5 au only the one-configuration result and con- 
elcerably easier to apply for any further caleulations. Wor R < 1.5 
the CE2 resonance mixes with the stabilizing configurations, so the 
entire eigenvector is presented for the stable Che root at each R in 
Table 7. The numbers to the left of the table indicate the space-spin 


configuration, with eleven spatial configurations and the two spin 





TABLE 7 
EIGENVECTORS FOR STABILIZATION OF on CE2 ROOTS 


Basis "OrbilGaliseror sCelwicnralions 


Configuration Orbitals 
l lo lo 60 
on Sues 
C how. So 
mec ea 
3 oOo 
| g uu 
4 lo lo 30 
ap wi Obl 
5 lo 20 60 
g uu 
6 lo 20 50 
ee aii aibl 
7 iio)” lemon 
S see 
8 Wonton: Sie) 
ee ae 
9 Uo Ml ephe C2) 
& & & 
LO To lo e360 
E& & be 
iE lonlosba 
& & & 
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Peublets. The odd indices correspond to the combination of a par- 
ticular spatial configuration and the triplet spin parent V/V6 
(2008-apa-Ban), while the even indices correspond to the combinations 
with the singlet spin parent 1/Y2 (afa-fam). The orbitals from 
Tables 4-6 which make up the eleven configurations are denoted in 
iaiple /. The eigenvectors con ane Pangest., cClhevenm Gomimarurar.on, ; 
stabilizing calculation are given rather than a smaller calculation 
foe the same stable root in order to show the relative mixing of the 
different stabilizing configurations. 

These results, and the results of the upper CH] calculations 
maemopresented graphically in Fig. l. The solid lines correspond te 
resonances of H,, and the dashed lines to Bee eetnon model calcu- 


2 


(ataons of the parent H, States. “NGveyliateds sa reclir On Meine a 


Lo M5 orbital as one of the major elements of the basis set in order 

immoeu rapid Stabilization of the core-excited resonances, the simsre 
Senripuration calculations for the excited i, 
Mieweonl to Che ecrossingwpoint 1s determined relative to tae 2e H 


States are Quiue food. 
fa 
eurve calculated with the same basis. 
, Cyt ; eo 
iets aiieo tO be tnereamunat the L revues ive CUrvevana tne L 


3 


eurve of e°- II do become nearly tangent at small R values. In fact, 
inom the Slopes in Fig. 2, it is possible that they may intersect. 
mhierefore, at two small R values, 1.05 and 1.2 au, the ad a9." and 
21 a configurations were mixed again as a check. The off-diagonal 
emements were ies smaller than the diagonal elements in both cases, 


Midiceating Liat wiepner Lhe non-crossing rile holds or not. the maxi 


of the two states is so small that if the nuclei possess any Kinetic 
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energy at all, they will follow the diabatic curve rather than the 
adiabatic, as discussed in the preceding section. 

For the stabilization of the two SP resonances of symmetry 
oN an entirely different procedure was used. Since the SP resonances 
eee physically considered to be a o, Clectron iim the tela or ei ke 
Papently polarized ground state, the first step was a full three- 
fenripuration variation treatment of the ground state of Hy The 
a and cat are thus not the same as those listed in Table 4, but 
Ghose which make the best ground state He One cern calculavioiee o IE: 
and a. is were cicseneds Wene r and x! as correlating Orbitals wired 
Minimized the ground state energy. This had to be done at each R 
Peat, in general, the orbitals with eee @ are much wniere 
concentrated about the nuclei than the orbitals of Tables 4-6. 

The first three-electron wave function was {6,¢,180 see 
oulsg 1s’ _ C4m2p,,x°2p/, } x jor1s + cis’ and the roots obtained 
with this calculation were checked for stabilization by adding the 

oF orbitals from Table 5 one at a time to the H ground state function 
mie opie first bracket. The trial funct@on above mixed Strongly with 
mie first three of these stabilizing confipurations, and tie roecs 
Sol began to settle down after the fourth excited oi orbital had been 
eaded, 

This was a very extensive calculation involving a maximum of 
twelve one-electron orbitals and as many as twenty-one spatial con- 
figurations, which with the two spin doublets made a 42 x 42 secular 


equation. For this reason, only the eneérsies and qualitativesfeauunes 


Glsune wave tunccions have been presenved. “Ine parameters oi all” the 





ihe 


@rpitals, the eigenvectors, etc. will be furnished upon request. 

As in the CE2 cases, addition of stabilizing configurations 
tended to drive the root corresponding to the SP resonance configu- 
ration up, inserting spurious and unstanle roots between it and the 
ground state and improving the ground state root. In contrast, CEL 
resonances calculated by the quasivariation method, when stabilizing 
femrreurations were added, tended to have their energies slightly 
Mmerered, if any energy change occurred at all. 

ii te Tiberestane vO Nove Liay anoluner rool in ae os sequence 
meeame approximately stable, at energies about 175 - 2.0 eV above the 
meston Of the four states Pai and sa discussed prev 1ousis)eebows 
Probably corresponds to a poor calculation On seme on (0° ED) 
mentioned above as a possible core-excited resonance. The singlet spin 
Mearent was definitely more important in’ the configuration, so it is 
also possible that the root could correspond to e Bot, however this 
fPeeelt to be léss likely osince The stabilizing cont tetra onesmend 
iemoesi Lhe rool up, not down. Further study 1s) necessary LO recclve 
mals DOInL, Since no configurations ovher sean E0024) - nd were 
miermuded in the trial function, and these canner be expected toncine 
aeecood calculation forvayCh state, even if Gestaple root Shoutdroceur. 
At present, there is no way of knowing whether the upper root would 
even remain stabilized upon addition of core-excited configurations, 
as it did not seem to stabilize as rapidly as the SP root. The semi- 
stable roots are given in column (f) of Table 3, and the parentheses 


are used to indicate that they are not considered final or completely 
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Pe2pilazed, buited® indicate the possible presence of a CE type 
resonance near that energy region. 

Finally, it should be mentioned that calculations of resonance 
energies, particularly by the quasivariation method, are fraught with 
unforeseen and hidden difficulties not encountered in stationary state 
miculations. First, in every resonance calculated, if the energy was 
plotted as a function of the "screening parameter" 5 of the outermost 
added orbital, pseudo-minima were encountered, a situation which had 
mo, occurred in the case of stationary state calculations with the 
diatomic molecule program used here. This is probably vecause 
variation of 8 shifts the maximum electron density of the orbital so 
that it coincides with the maximum dene Of higher orbitals wieven 
imfjouch the orbital exponents n and m aré fixed and do nov represent 
the nodal properties of the higher orbital properly. Hach coincidence 
of electron density maxima would then produce a pseudo-minimum. An 
eeomple Of this is the wguasivariatioen calcullatlonge: suhe sec 1 
feesonance. At R= 3.O0-aus~essentially the same enere, can ero vaimed 
Perch 2. (.03056,.0;12-11) or 9..(.0632,0;42-11) and fixed #, = lo, and 


3 


B, = 1m - Both . OorDipeliswexhibit false@manima Tor Miweer 6 values, 
The second difficulty is with convergence of the two-electron 
integrals (ij)d/rik2) when very small & values are used, particularly 
mi the higher orbitals, Within the limitations set by the program, 
firs problem 1s purely numerical, caused by the loss of Signaticant 
figures upon subtraction of large numbers and exponent overflow, 


underflow and wraparound, and can be surmounted by the use of scaling 


and/or double precision arithmetic. 
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itd. tiiere 1S) the anount Of paysical and chemical» antuition 
required to obtain good results. A poor guess of the major resonance 
Boniigurations, or a poor trial function for them, may either greatly 
increase the amount of work necessary to get agreement with experiment, 
or yield completely meaningless results. The ability to calculate 
widths and energy shifts would greatly aid in the problem. 

The consequence of these difficulties is that each resonance 
calculation must be approached as a completely separate problem 
meaquiring special, individual handling. There can be no doubt of the 
fact that the quasivariation and stabilization methods work, as the 
aereement with experiment reported here shows, but for these and any 
@ener calculations of resonant states, pours Caui tone oe tece se cia, 
With a well determined basis set, and based on our experience with 
ig@tm methods the stabilization method appears toybe preferable for 
reasons of both flexibility and safety. 

The calculations were performed on the Honeywell 800 computer 
of the University of Southern California Computer Science Laboratory, 
and the LBM 7094/7044 DCOS equipment of the Jet Propulsion Laboratory. 
ine authors would like to express their appreciation to both 
institutions for making available the large amounts of computer time 


which were necessary for these calculations. 
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CHAPTER IV 
SUMMARY 


Quasi-stationary methods of calculating the energy of core- 
Peeived and single particle resonances and autolonizing states of 
atoms and molecules are presented both from an intuitive, physical 
approach and from a formal, rigorous standpoint. 

The methods are applied to resonant states of H» and in the 
Pmecess, & general computer program for calculating the energies and 
wave functions of diatomic molecules which incorporates DOUN eras 
normal variation method of quantum chemistry for bound states and 


gquasi-stationary methods for resonances has been developed. 


The results of the H- 


5 calculations are compared to electron 


EpeauLvering Experiments. 

The diatomic molecule program developed is ea modernizatieneor 
mmevious methods including Several new options, and “extended so mthar 
colemlalions can be “dene for “up to twenty electrons] tv ais ful s— 
eussed in a set of Appendices, and a sample calculation on NO* is 
presented as an illustration of the applicability of the program to 


feateulavions on large diatomic Systems. 
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APPENDIX I. MATHEMATICAL DETAILS 
Wave Function and Coordinate System 


The diatomic molecule program discussed in these appendices is 
a generalized program for calculating the energy of a diatomic mole- 
cule within the Born-Oppenheimer approximation, using the variation 
method with optimization of both linear and non-linear parameters to 
in prove the energy. The mathematical procedures used are an extension 
and generalization of the programs of Harris (1960), and Taylor and 
Harris (1963, etc.), and extensive use will be made of these references. 
The inclusion of the full mathematical treatment, rather than only 
the outline of its departures from the above methods, is for the 
purpose of completeness, so that these appendices will contain all the 
information necessary to understand and to use the diatomic program. 

The total wave function is composed of linear combinations of 
products of one-electron molecular orbitals. These orbitals, B., are 
expressed in prolate spheroidal coordinates, therefore the discussion 
will begin with the coordinate system and definition of the B.» and 
proceed stepwise to the definition of the total wave function. 

The prolate spheroidal coordinate system is based on two foci, 


A and B, where the two nuclei are located. The foci are separated by 


the internuclear distance R. The coordinates of a point Pj are 





Ell 
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?, = azimuthal angle about the internuclear axis. 


— 
ae 
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The quantity Os» uhe distance of Py irom the inGernue lear wants, 
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The volume element and some differential forms in this system 
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The one electron spatial function B. is chosen in terms of these 


Peordinaves as 


(Ge) 


F n> mj Vs 
exp |- 5.6. - an, tivi@s|- 6, a . (20,/R)! 7 (6) 


O, (5,50, 3n mv.) 


where 5. , are adjustable non-linear parameters which will be 
determined by the variation method. The 5 parameter is a measure of 
the extent of an ellipse of electron charge density about the two 
nuclei, and must be positive to assure that the function goes to zero 
as & goes to infinity. The ratio a/& gives a measure of polarization 
along the internuclear axis since @ = 6 centers the electronic charge 
on the nucleus at A and @ = -6 centers the charge on the nucleus at B. 
Intermediate Q@ values allow improvement of the molecular orbital by 
allowing unsymmetrical, polarized charge distributions. 

The exponents nN,» M,, are restricted to positive integers or 
zero, and are used to form higher orbitals. Non-zero Mis introduce a 
node in the perpendicular bisector plane of the internuclear axis. ve 
is the azimuthal quantum number, and can assume any integer value, 
including zero. The (20,/R)! 3 factor has been included to provide 
regularity on the internuclear axis, and because it provides simpli- 
fying cancellations in some of the integral expressions to be 
introduced later. 

Atomic orbitals, in particular Slater type orbitals, can be 
expressed as linear combinations of the elliptical functions B. with- 


out considerable difficulty, elthough the number of terms necessary in 





8O 


the combination goes up quite rapidly as n increases. For example 


a a(8, 
I, 
es. = 2 
Teal 

es, = aR 


SREO IL TS = 8(5,-83;000) 


R [9(5,8;100) + 8(8,8;010) J (7) 


[2(8,-8;100) - 8(8,-8;010) ] 


The remainder of the series can be generated using the fact that 


from (1) one can express 


In the homonuclear 
ungerade symmetry of the 
There are two methods of 
the total molecular wave 


taking the proper linear 


AL 
ea equa) 
(a) 
Y= A R(E— - ) 
b 2 q 


case, the problem of obtaining gerade or 
electronic wave function also presents itself. 
accomplishing this. One is the passing of 
function through the center of inversion, 


combination of the function and ats iiverted 


counterpart. This method will be discussed later, after introduction 


of the total wave function. The other is to symmetrize the individual 


molecular orbitals B.> actually defining them as oe or eC etc. 


The symmetrization 


of the molecular orbital Bb. 1s accomplished 


by inverting it, and taking @ + ig with the plus sign for gerade and 


—>> 
the minus sign for the ungerade case. The i represents the inversion 


operation, which is equivalent to a reflection in the perpendicular 


bisector plane of the internuclear axis (interchange of A and B), 


followed by a rotation of pi radians (p> + n). 








owl 


Examination of the function @ (Eq. 6) shows that this inversion 


is accomplished by replacing yn by -yH and m by » + x, and therefore 


BU OC 


g = (8 ,asmnv) +: (ep) 6(5,-asnnv) } yo (8) 


where k = p +m + ae p representing the parity g or u desired, with 
feeven for g and odd for u. 
Since only certain portions of the function are affected by this 


operation, it 1S convenient to rewrite the symmetrized orbitals as 


g.®"(6,a;nmy) = et E™(op/r)!¥| ie | ik. (-1) "ein 1 
(6") 
to eliminate exceSSsive computational labor in the actual calculation 
of the spatial matrix elements. 
Now let us define a primitive function o.(z) 


Mt > 
ME ei) (9) 
=i 


Es 
?.(r) = 

where N is the number of electrons in the system, ca is the position 
vector of the es elleeuron, r symbolizes collectively the coordinates 
of all N electrons. The index j will be used to denote a particular 
spatial configuration after configuration interaction has been 
introduced. 

A single spatial configuration is the antisymmetrized primitive 
function. Combining it with a spin function @, one obtains the 


expression for a single configuration 
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where nt is the a eigenfunction of = and G represents the spin 
coordinates. Usually in this program one spatial configuration is 
combined with several spin eigenfunctions comprising a complete spin 
basis, but depending on the number of closed shell electrons selected 
in the problem. A more limited spin basis may be used, or a specific 
Spin eigenfunction may be selected which most corresponds to the : 
physical situation. Once a specific problem has been selected (i.e., 


lowest lying 3 


= state of the molecule AB), the number of spin eigen- 
functions will be fixed, so for Simplicity the word configuration will 
be used to refer to a spatial configuration, realizing that this 
spatial configuration may be combined with one or several spin 
munetions. | 

Defining d as the number of spin functions chosen for a par- 


ticular problem and m as the number of configurations, the total trial 


wave function for the ie state of the system becomes 
> > = L idee 
¥, (59) = me » Che Slhcoee a) 0) (Gia, 


where the Ci are the components of the ie column of the eigenvector 
matrix determined by diagonalization of the secular equation. 

It 1s now possible to show the method of producing g or u 
Symmetry by inversion of the entire molecule through the inversion 
center. 

The linear combination 


? 
e,¥ + cai ¥ (12) 
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eam be considered 4s a two configuration problem, and the total 


secular matrix will be 


(yy)  (w*HY) 


ls 
(iv*Hy) (iv*diy) el 


which is a supermatrix of m dxd blocks, and is of the form 


A B 
( (14) 
Bt UA 


since the matrix is Hermitian and (iv*H;yv) = (wtHy). 

If this is diagonalized for the energy, both g and u solutions 
will be obtained due to the form of (12), but it is well known that a 
matrix of tne Subnet (14) Can be diagonalized by tne unitary, 


transformation 


iE 
We i>) 
: -1 1 
and applying this to the supermatrix (13) we obtain the block diagonal 
matrix 
(fy + iv]*H [vy + iv]) O 
iH 
5 a) 
O ((y - iv]*H (vy - iv]) 


or, taking advantage of (14) again, its simplified equivalent 
(v*H (y + i¥)) (16) 


Thus we see that Eq. (16) is equivalent to actually taking 


a = 42/2(y¥ + i¥), but the amount of labor for calculation of 
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(16) with the proper sign is only about 3/4 that of calculating and 
@iaponalizing either the g or u block of (13"). 

The program itself calculates the blocks A and B and then takes 
the sum or difference as necessary to select the desired symmetry 
for the molecule. 

iierie verse Of 2 States, if pl oOrvlvtals are samelnded piuc- 
minus symmetry must also be considered. This can be handled by 
replacing va by “Vv in each pi orbital in a second configuration and 
repeating the calculation, and a procedure completely analogous to 
the above for selecting g or u symmetry can be used to select only 
the plus or minus linear combinations. Also, one could define a By 


by replacing the exp(ivm) factor in @ by 


[exp(ivp) + exp(-ive)] (17) 


mi the pi orbitals of a configuration. 
Spin and Antisymmetry 


In the original diatomic methods (Harris, 1960; Taylor and 
Harris, 1963), spin and antisymmetry were handled by the method of 
representation matrices presented by Kotani, et al. (1963). This is 
an elegant group-theoretical approach in which a complete set of spin 
eigenfunctions is used as a basis to construct irreducible repre- 
Sentation matrices, suitably antisymmetrized, for the ake order 
permutation group O.: The method is completely general, and is 
particularly adaptable to the case of non-orthonormal spatial orbitals 


in the complete open shell (different orbitals for different spins) 
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treatment of the diatomic program. 

Unfortunately, the number of mathematical operations in the com- 
plete open-shell treatment increases astroncomically for larger systems, 
due to the n! permutations in the group o and the dimension of each 
of the set of n! representation matrices. For modern computers, the 
practical limit appears to be between six and eight electrons. In 
order to study larger systems, some sort of approximations must be 
made just to reduce these numbers to something in the realm of 
manageability. 

One possibility is the closing of some or all of the shells 
(double occupancy). This has the effect of reducing the dimension of 
the representation matrices to only the invariant part corresponding 
to the remaining unpaired electrons as is outlined in Kotani, and of 
reducing the number of permutations necessary, since double occupancy 
causes Linear dependencies in the group o: If there are r paired 
orbitals, the number of permutations is reduced to nt/2". As an 
interim measure, this is helpful, but as n continues to increase, even 
closing of all the shells fails to reduce the problem to one of 
manageable size. For example, 161/28 is still of the order 14!, and 
even though the representation matrices are of dimension 1 X 1, it 
is clear that some other approach is necessary. 

Another alternative is to introduce spatial orthogonalities, 
and by knowledge of where these orthogonalities exist, eliminate 
those permutations which can make no contribution to the total secular 
equation. This procedure can be used independently of the shell- 


closing procedure, or they can be used in conjunction. In the 
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16 electron case mentioned above, by orthogonalizing ten inner shell 
@rbitals in closed pairs, and leaving six open-shell electrons for 
bonding, the size of the problem is reduced to somewhere between the 
seven and eight electron complete open-shell problems. This is still 
rather large, but well within the capabilities of a modern large 
scientific computer. 

Thus one of the aims of this research has been to revise the 
Spin-antisymmetry routines of the older versions of the diatomic 
program to include options which allow a range between complete open- 
shell and complete closed shell spin treatments, using any ortho- 
gonalities which are introduced or exist naturally to reduce the 
problem. 

Since the formal treatment is fully presented in Kotani (1963), 
the remainder of this Section will be mainly concerned with the 
actual mechanics of applying the method with all the options as it is 
done in the edition of the diatomic program presented here. The 
notation used is chosen to be as consistent with both Kotani and 
Taylor-Harris as possible. Unless otherwise noted, arguments of 
functions have been Suppressed, and the coordinate number is indicated 
ie position (i.e; 6(r,) pe) 625) = 90,9). 

Let the two eigenfunctions of the one-electron operator s, 
be denoted as O and f, with eigenvalues 1/2 and -1/2, respectively. 
Then eigenfunctions of the many-electron operator 5 can be represented 


Sepoimole produets ol GQ and 6, 


Q=AQ+++ BOB: . (18) 
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With eigenvalue M. In general, there are 


a (19) 


menearly independent simple product functions 2 for a given M value 


in an N electron problem; and the desired relation 


N/e2 N 
iL 
M=-N/2 5 N+M 


Ties funceioms are nol generally eleenrunceions or ce which, 
ike 5. commutes with the Hamiltonian, but a simultaneous eigen- 


mune tTLON of s and S. ean, be “expressed “as 


N 
Os = Ye c. 0, (20) 
a 


N N 
Genes Ss) =" en Ss 2 etncmaons 


t] 


There are pl 


S 
N N : @ 
q k = 1,2,°°:F, for given eigenvalues S(S+1) of S and Mof S.. 
5,M,k S =_ —z 
In the diatomic program, the coefficients in the above expansion of 
om MR ere Chosegeeneratedsby the method of senecailesg calseonstrucuion 
pee | 


(Kotani, 1963) which uses the vector coupling formula and the functions 


eeand & FO produce 


gl 7 S-M+1 ,N-1 of Sith QN-1 : 
S,M,k eS +2 $+$,M-53k eS+2 “$+5,M+S;k 











88 


N S+M N-1 S-M .N-1 N-1 
6 = ~——— § al & ———— § “Rp |e + h 
Suk ~ YS “saudi YS “s-z,medsn’! Sid 
h = ee ody 
S73 
(il) 


iiss nova unique choice of oN but any other choice may 


SM, k? 
be expressed as a unitary transformation of this set. Katsura (1963), 


N 


has published a set of all Og M, k 


Up to Sixeelecctrons, Vand Welt am 
(1966) has developed a rapid and efficient scheme for generating all 
Bx which is limited only by available storage space in the 
computer. 

For the case of d and m both equal to one (one configuration and 
one spin function), so that the total wave function is represented by 


Eq. (10). The matrix element of an operator 0, which acts only on 


Spavial coordinates, is 
(A, [*9; J JO] A, [*,°,] ) (22) 


and using the Hermitian property of the antisymmetrizer plus the 


; e 
relation A” = yin A, = 2 (-1)P p, (22) becomes 


Dy (-1)* (6,6, || P|, 6,|> eam 


p 


where p is the even or odd parity of the permutation P. Since 


P= PP, Eq. (22') can be separated to obtain 


> (0(7)[olo(P,2)) (6(8)| (-1)?6(P_3)) (23) 
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where the (ais) has been included in the spin integral for convenience. 
Calculation of the spin integrals yields a set of coefficients U(P), 


so that the matrix element can be expressed as 


22, 2){e1 9 P0) 2 u(P) o(P) (2h) 


with the definitions 


If 


(o|0| Po) 


I 


O(P) 
(25) 


U_(@) (-1)? (e|Pe) 


where the permutation is on either the functions or the coordinates in 
both integrals. 

Additional spin functions may be introduced in which case the set 
U (P) becomes a matrix of dimension d, where d is the number of spin 
elgenfunctions in the basis. Thus for one configuration, the matrix 


element of the operator O becomes a dxd matrix of the form 


> U(P) O(P) (26) 


P 


where 
u(r) = [u,,()] = (-2)? €, [P2,) (27) 


The secular equation for the energy is then 
2 U(P) [H(P) - ES(P)] = 0 (29) 


If the complete set of spin eigenfunctions is included (d = Fo) 


the matrices U(P) become unitary, irreducible representations of the 


permutation group, and are exactly the U(P) defined in Kotani (1963) 


and the Taylor-Harris papers (1960, 1963, etc.). 





When spatial configuration interaction is included, the 
definition of O(P) can be expanded to 0. (P) = (,[0/Po,), and the 
total secular matrix (29) becomes an mxd matrix, composed of dxd 
blocks for each Hy (P), S . (P)3 jie (been cia 

Given the sets me it is now possible to generate the repre- 
sentation matrices U(P) systematically, thus completely separating tke 
spin and space parts of the problem, since once generated, the U(P) 


are general for any problem with the same N,S, and M. From Eg. (27) 


Ue) =6GE)\- Ce, (ie 


and therefore for an arbitrary permutation P, the associated U(P) may 


be generated from a single rectangular matrix C, which has Nn, rows and 


d columns stemming from the expansion of a. 


SM in terms of the 2. 


Wea. 20). 


Thus the p,v elements 


d d 
Uy fP) = Ped C59, Ip > c, 8, ) 


Jal 


d d 
(-1)* 2 2° ip Cp yy <9; 195 


where se = Opa? which is a member of the complete 2 set, and due to 


the orthogonality of the 2 
a) a> 2 canes Pv SPs 


iis 1s Yepeaved for alli, 5) lo ed sv ed colmake the full 


matrix. The a's and B's are easily represented in a digital computer 


Sy 
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as O's and 1's, making the operation He and the checking of or P3 
3 
relatively simple. 


Poredetrenmcaved OG. Dacis, Gepending® On Ge Tune Lions whiten are 


k 
omitted, certain rows in C may have only zero elements (Ci = (Ole 
1<p<d). If so, those rows are omitted, and the multiplication by 
zero handled implicitly by the Kronecker 6 function for computational 
@onvenience. 

Since it is now possible to form the matrix U(P) for any 
Pupinrary FP, the next requirement is a method of generating all the 
permutations P of the group o,. Harris (1960) devised a simple 
algorithm for systematically generating all N‘ permutations of the 
it group, pee mathematical operations and scan time. 

Consider the arbitrary permutation P. It is desired to generate 
mie succeeding permutation Q@ in as few operations as possible and such 
Gian @ 1S nob a repetition of some permutation generated prior to LF. 

Read the numbers in P from left to right until the first 
increase occurs (i.e., first number which is larger than the one to 
its left). The position of this number is called the pivot position. 
interchange the mumber in thespivot pesiticn with Chesiarces pagumoer 
to its left which is less than the number itself. Simply copy (leave 
unchanged) everything to the right of this position. Then arrange all 
Demers to the left of he pivot pesitiems in ccending yorden whe 
pecult 18 the permutataon.Q. Ef there wasano “inernesce wall pernu— 


tations have been generated. 
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For example, consider P represented by i/l j/2 k/3 L/ “S82, 
wen ft = |} < k > J, and i >k. The pivet pesition is 3, occupied by k. 
By the inequalities, j is the largest number less than k, so j and k 
mec interchanged and £....s are left as they were. Then i is greater 
than k, so they are rearranged in ascending order, and the resulting @ 
representation is k/1 i/2 j/3 41/4 ...-s/n. 

The parity is easily determined by inspection of the canonical 
number of the position immediately to the left of the pivot position, 
eee ft Cpiy-1 [Mod. 4] is less than two, the parity of Q is 
Opposite to that of P, and if it is > 2, the parity of Q is the same 


as that of P. In the example above c = 2, and 2[Mod. 4] = 2, there- 


Pel 


fore Q and P have the same parity. 
Table 8 gives a full numerical illustration of the results of 
the above techniques for the case of a complete open-shell three 


and 63 and six 


1 1 ) 
sayl 5,3 2 


permutations in the order they are generated. 


electron doublet, with two spin functions zr 


The above procedure is general for any number of 9 functions, 
but if closed shells are introduced and the dimension of the @ set 
reduced concomitantly, the U(P) become linearly dependent. (The 
expression "closed shells” is to be interpreted as strict double 
occupancy.) It is still possible to generate all the N! permutations 
and U(P), but a much more favorable computational situation would 
occur if a way of generating only the linearly independent subset of 


P and U(P) could be determined. 
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TABLE 8 


SPIN-PERMUTATION MATRICES FOR THE 
OPEN-SHELL THREE ELECTRON DOUBLET 


3 ua 

0 a ECOB-GEBO- 
Ba (2008 -aBa- Ba) 
3 iL 

0 = O- BAO! 





321 (-) 2( z 


For r closed shell pairs of electrons, the @ basis becomes 


i Gg 

N _ .N-2r ‘ Je 

ahi: a se M, k (aB-Ba) (@B-Ba) 7 (oB-Ba) ( 2D 
eee (30) 





on 


and there are Ni /2° linearly independent matrices and permutations in 
the subgroup of Oy: The generation of the U(P). is the same as before, 
except using the reduced @ basis so the problem becomes one of 
generation of only the linearly independent subset of permutations. 

The Harris algorithm can easily be adapted to this problem. If 
an auxiliary permutation vector R is introduced, such that the vector 
P contains the actual function (coordinate) numbers, and R contains a 
set of dummy indices in which the closed shells are represented by 
duplicated numbers, a linearly independent subset can be generated by 
doing the permutation-generating scan over R, performing the same 
operations on P, and using P to generate U(P) and O(P). 


For example, a six electron triplet might have the form 


Z 
ye (ap+Ba) (aB-Ba) (aB-Bar) 
and the corresponding initial P and R would be 


jee eee aie 


Rican ee ae 


As an even more complete example, consider the three electron doublet 
with one closed shell pair, whose spin eigenfunction is om > and 

RE 
whose P and U(P) can be represented by the Uso elements of the six 
permutation matrices of Table 8. Here it is seen that of the set of 
six permutations, if electrons one and two are the same, 123-213, 
132-231, and 312-321 are equivalent, and moreover, that the 2,2 


elements of their matrices are also identical. Thus Selection of three 


of the non-equivalent permutations with their respective matrices would 





SB. 


only have the effect of multiplying the total secular equation by a 
constant, in this case 1/2. 

Setting the initial P and R as 123 and 113 respectively, the 

sequence of permutations becomes 

eS We esi) Se) 

R- 113(+) 131(+) — 311(-) 
so we see that the search over R will produce 3!/2 = 3 non-equivalent 
permutations and using P to generate U(P) and O(P) will produce the 
desired secular equation. 

The preceding techniques can be applied automatically by the 
present modification of the diatomic program. The case of orthogonal 
functions is slightly less easy to generalize, and thus requires some 
planning on the part of the user before proceeding to the computation 
stage. No modification of the program is necessary, only careful 
preparation of the input data cards. 

A completely orthogonal spatial basis would have non-zero ele- 
ments only in the case of pairwise interchanges and the identity per- 
macation. 

For any other case, one must consider the spatial orthogonalities 
and the coordinate operators. Clearly, the set of permutations used 
must include all possibilities of one or two electron operators con- 
tributing to the secular equation. Thus, any permutation of a con- 
figuration which produces no more than two zero factors in S(P) must 
be included, and conversely, any permutations which may be seen to 


include three or more overlap factors which are zero may be excluded. 
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iia many cOnfiguration calculation, the configuration with the mini- 
mum number of orthogonalities controls the exclusion of permutations. 

This option has been included in the program by allowing the 
user to generate a Set of permutations particular to his problem by 
hand and instruct the permutation generating routine of the program 
of the specific limitations placed upon it. The program will generate 
any additional permutations as required and produce the reduced Set 
of U(P) delineated by the permutations. 

As an example, see the NOT calculation reported elsewhere in 
this work (page 140ff.) which was included specifically for the 
Meepose Of illustrating the use of this option. 

A word of caution is neccesary here. When a complete spin basis 
is used, and the full N: set of matrices U(P) is thus generated, the 
ordering of the one-electron spatial functions in each configuration 
is energetically immaterial, since re-ordering is a unitary trans- 
formation, and the eigenvalues will not be changed. However, the 
elgenvectors will reflect the re-ordering transformation, and they will 
be changed accordingly. 

On the other hand, if a truncated spin basis is used, whether all 
N? U(P) are generated or not, the eigenvalues will be dependent on the 
order, and for best results, the spatial functions must be ordered to 
correspond physically to the ordering of the @ and B in the spin 
functions. 

If one takes advantage of the orthogonality options, obviously 


the functions must be ordered to correspond with the set of hand- 
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reenerated permutations, and the ordering must be preserved in all con- 


figurations of a multi-configuration problem. 
Matrix Elements 


To produce the secular equation (29), the spatial integrals H(P) 


and S(P) must be evaluated. From Eq. (25), these are recognized as 





85(P) > (*51"*,) (31) 
and 
H, (BP) = (:,§,) (32) 
where 
Ve n n 
A'B 1 
SS ee > ae (33) 
1=L 1 Sl) 
and 
Z, Z 
oe A B 
h. =-=V, -—-— (33') 
, an iA “aR 


can and a are the charges of nuclei A and B, Separated by the 
distance R. The distances ee are the distances between particles x 
and y. Atomic units are used, with energies in Hartrees (1 au = 
27.21 eV), and length in Bohrs (1 au = .52917 A). 

Since no term of H# includes the coordinates of more than two 
electrons, H(P) is reducible to sums of products of one- and two- 


electron integrals, and obviously S(P) can be expressed as a simple 


product of one-electron overlap integrals. 
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Thus the matrices 


(ilhla) = fo,*Gr,) m9, (F dar, (3H) 
Gin fore) a(t )ar, (35) 


and 
as > > Ae > > > > 
Grae ice) - foe) 03°.) Es g(r, 8, (r,) ar, dr, (36) 


can be calculated and stored, and the necessary elements combined to 
make H(P) and S(P) for a given P. 


Expanding in terms of the coordinates (Eq. 1), the one-electron 


E 


integral (ilh|j) separates into a kinetic energy integral and two 
nuclear attraction integrals, 
aay o 


~=v* )-§ -2p(! ) 
: | a (E°-n") 


2 (37) 


With the form of the wave function (6), all the one-electron 





(| hee (: 








+e Gio Be) (: oe 
Grea) 





integrals may be expressed in closed form in terms of the functions 


(Kotani, 1963) 


A (8) : t exp[-8t] at (38) 
1 


ug 
B(@) -f t exp{-ot] at. (39) 
-| 
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The integrais for non-zero v may be compactly handled if the defi- 


nitions of (38) and (39) are extended to 


Then, by direct substitution of the function (6), followed by 


(38°) 


9") 


separation of variables and integration over 9, it may be found that 


3 
(15) = SE o(vz,v,) [ay (6) BY(a@) - a¥) 














WMiere N =n. t+n., =m. +m., 6 = & + 6, GQ= a. 
at J a J a j 1 


The kinetic energy integral is more tedious, but after some 


manipulation can be converted to the form 


+ fop\l" 
(31- = v3) = - - suv) f a (=) 


Vi s-ok oe ea or eee (2) 
i i j i i j 0 


Ny My 


3 
} | R v x 
(: (¢2_42) ' = “To 8049¥5) Aner 8) B,C) 
; Se oe ; _ AR? S(v. v.) A‘ (8) B 
(4°17) J My ule J Yn 


Ss 


(0)] 


(a) 


. and 


5,°0,9| 
aie 


(40) 


(41) 


C2) 


where o.- = ¢ 7 exp[- (5,6 ta.) ] , and after substitution of (6) 





LOO 


becomes 


(i[-5 9°15) 


3 
i a 


D 


°) AY (8) } BY (8) (43) 


V / 
+ 28/n’ A768) + (n-n 


Vv 


eV 
: fe B(a)-2(a/m xe) se 


i (c7) #200/m/BY 5 (1) +(m-m’°)BY_ (a) a” (8) 


TO , 2, 12 a -l lt -1 
fr tn-8 oar ater’? haY (6) BY (a) ty yeni3(8)B (a)-Ay cot, (0)t] 


eee 1. = hoy mM’ =m. = m., 6 = 6. = 60, and a” =o weno 
L J 1 J 1 J 1 J 
fmiescoectficients of A (8) and B’ (a) vanish when negative indices 
would be required, and correct results are obtained when these unde- 
Vv 
fined values of AL and Bo are assumed to be Zero. 


The symmetrized orbitals (6') are easily included by noting that 
ae 
een eu) 
: ky+k >. . ki. : Ki;. : 
mee) ly FC oa) (1) ots Gs ecu r eles 
where i’ is the function @(5,-o;nmv). It is clear that 


BY(-0) = (-1)"8Y (a) 


therefore it is possible by collecting terms to define the auxiliary 


integral 
/y g Vv k-_v 
B, (a,07) = 28[0,m4ke, +k, (mod »| {8¥(0) + (-1) JB (a7) \ 
end since in each of the integrals (40-43) all terms involving m are 


of the same parity, the Kronecker delta function can be factored out, 


leaving the integrals in exactly the same form, except using 
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B’(a,0’) rather than B (a). This method is not unique, but it has 
the advantages that it requires less memory storage, and reduces the 
amount of actual computation by avout one-third. 


The A (8) are obtained by the recursion formula 
= e 
Be) (Se (e) (4b ) 


with A (8) = ge e >. The B (a) could also be generated recursively, 


using the formula (Kotani, 1963) 

B (a) = -  (-1)" oiene ee nB_, (a) 
but Corbato (1956) has developed a method which does not involve sub- 
imaction Of terms OL Similar magnitude and therefore is more accirace. 
The B are calculated in Conjunction with réelaved Tunctions vseqean, 
mae two-electron integrals. 

Because all the orbitals in this problem are of a single para- 
ferric foum, there is bubeonertypevor c-eleetron inteanal in cur verk 
meyer than the usual. plethora of hybrid, coulomb, exchange, sandgoac- 
Sewer integrals. To evaluate the 2-electron integral Gayl seis) we 


use Neumann's expansion (1878) of l/r 





ee 
00 Li 
a 
Re ie O Che hole 
or 5 » ) er ee ee ese] 
1 
u=O O=-4 


(45) 


x Beja tee lg yee! m,yexo[ iol, -9,)], 


mere Ey and tlemean che larger, and Smaller, reapeceively. ot 
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e1 and Eo° p | and a, | are associated Legendre functions. 
On substituting Eq. (45), and introducing the wave functions, 


Eq. (6), we perform first the integrations over @, and @,. If we let 


a | ee 
K i ees 


Vy 294 =o (46) 


the m dependence of the integral is exp i(o, +0), + i(9,,-9)®, cet 


integral vanishes unless 
Se ap Ol Oh (47a) 


in that case the sum over co collapses to the single value (or its 


negative) 
o= Jo.) = Jopl (WT) 


Henceforth we shall use o to mean the value given by Eqs. (47a)-(47b). 
The m integrations then cause the removal of the o summation, and the 
change of the lower Limit of the w summation to o, while contributing 
Pomc! ri/kf) a factor bn 8(0, ,-0,). 

We next consider the factors involving Cea and Cone When 


the contributions from all four wave functions are collected, we have 


(..°- )(.-n,") 2lviltalv«| (..a) 0-5) Liva| +l ve| 


When we do the € and j integrations, we will find it convenient to 


Write this as 


(a : ( n°) ho" (. :) ( ) — 
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since the 1/20 power can be naturally associated with the Legendre 
functions. In general, the t value to be associated with any two 
orbitals, say 8, and B will be Tt = 1/2(|v f+lvt-[v-v 1), and in 


marvicular 


1 =3 (Iv sll ygl-2) 
1, = 5 (Iv 111-2) | (48) 


Note that these definitions imply that t is a positive integer or zero. 
We are now prepared to consider the integrations over G1 and Boe 
It is evident that the integrand cannot be factored to separate the 
coordinates (E55 nm) from (ee No) because of the presence of the 
Legendre functions involving €+ and €_. This situation can be 
Bifilcviated somewhat by using a slight generalization of a formula 
published by Ruedenberg (1951). The formula we need can be expressed 


as 


(-1)°(2u +1) fa, fh 240,229,666, 
ik it 


f axg, (x)6, (x) (49) 
i 


where ge,’ (6) and a.” (x), y = 1 and 2, are related by 
a6 
G” (x) = (2H +1) [las 3)!/(u-0) J 


xX 
x Pp 209 0 - ial f agp (E)e7(E). (50) 


a 
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This transformation may be applied whenever all the operations 
indicated are meaningful. The exponential dependence of the wave 
functions is sufficient to guarantee the utility of Eqs. (49) and 
(60). 

Summing both sides of Eq. (49) from p = 0 to 0, we have, with 
suitable definition of e,’ (E), an expression for (ay cee Exam- 
ining the result of the substitution of Eqs. (45) and (6) into 


(ig{1/r{kf£), it may be seen that 


1 
@” (6) = Eou-0):/(uso):} ff an(E--n- Jexp(-8 8-0, 1) 
an 


x PL°(n) (B21) SOUP) #87 (1g? (0/2) #7 gy ty (51) 


: i _ re ; ee 
Here, rT o. Os 5, = ve with corresponding definitions tor 


a, 1, and nl, . The n integrations appearing in Eq. (51) are all of 


the same general form. Define 


i (m,a,7) = 4(2u oi) [Cu - o)t/(u + a) 
1 
«f anP (1) exp(-an)n (1 - yer (52) 
Sil 


and Eq. (51) may be rewritten 


(2u41) [(ut0)!/(uto) Tbe 7 (€) 
m m m i : BS 
n 2 GJ 2yaT 2) b 16) 
= exp(-6, 6 )E Y(E--1) re i (m,,@,,T,)-i (n, +2,0, ,,)| 


The integrals defined in Eq. (52) are related to spherical Bessel 
functions, and also to the B (a). The B_ can be expressed (Corbato, 


1956) in terms of the i as 
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‘s 1 
/ ray 2 ea eg 
> ee (54) 
Gia) ee Gane 
H=0,1 


where o = 0, or if v # O, since o is necessarily zero for non-zero 
one-electron integrals, a Ye = T, so the binomial summation of 
(39') can be avoided by using the relation 


y ee zs 2 
B (a) = >) Ht i (0,a,7) . (54") 


CoAT atte IY 
n=0,1 


In both expressions the prime on the summation sign indicates 
taking only every other term such that the highest pp value =n, and 
the double factorial is n!! = n(n-2)(n-4)...(2 or 1). 

For the symmetrized one electron orbitals, the i” (m,0,07,7) 
become 


oO Kae FO 
i (m,a, a’,t) = {5,,°(m,a,1)+(-1)54,,%(m,0”,1)} 


“Oro moda (m+kytk; +1-0] 
(55) 


Since 


O iid eG 
P = = 7 P 
en ee) 1, n) 
The next step is to substitute Eq. (53) into Eq. (50), and per- 


form the — integration. When doing so, it 1s convenient to define 


another set of integrals 


e Oo eae a bo 
1 (n, 657%) a IP. (x) (x <p) | (56) 


. ‘ 
xf abp,°(¢ Jexp(-Be )E™(E°-1) (/7)** 
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in terms of which Eqs. (53) and (50) yield 


94 0 5 20 
G X) = K - ego 
tt ( ) UL (n, , yr ty Xi, (m, ,@,,,,) 


Oo Pai 
- K (ny 8,7, x) 4 (m, + 2,07, ) : re. 


To show more clearly the role of the wave functions in Gg” (x), we 
shall from now on use the notation (u} G(x) |v) to mean a form like 
y F 
Ga (x), but with ny replaced by n +n, etc. 
Recapitulating the earlier discussion in more explicit form, we 
obtain the result that the two-electron integral has been reduced to 


the form 


3 Z CO oO 
(ag|2/elk2) = (= (Joe ox) (316,00) 14) (316,00) 12). 
Jk =O 
(58) 


The integral over X in Eq. (58) is now to be evaluated numerically. 


Using values of the integrand at points X 3 the integral can be 


written 
feng |L/2| Kz ) 00 
a (2) (Rove, odor, 2d (2/6 (x, )]k) (GIG (x, ) 14), 


where a is a weight factor. Defining 
L 
ree re eee : 
(ww? Gla (15) = GIG 13), (59) 
Eq. (58) becomes 


Gjli/r|ke) 


3\¢ = 
; (=) @ Teme Z. 2 GIG, POGIG TH) = 
K  {4=0 
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The number of terms one must keep in the sum over wu may be 
estimated by studying the behavior of the K and i integrals, Eqs. (52) 
and GeOR By far the most important quantities affecting the con- 
vergence of the sum are found to be |x, | and la], Pano cmee vale see 
which correspond to slower convergence. In fact, for @ = 0 the series 
reduces to a finite number of terms. Since the speed of computation 


depends on the number of terms, we actually compute a number chosen 


for each integral after examining its |o%, | and |, |. As a rough 
iadication, we use an average of 13 terms for | or] = 10, and 19 terms 
for the maximum value, jal = 19.9. 


Details of the calculation of the integrals i (0,7) and 
K  (n,8, 7%) can be found in Appendices II and III of the Harris 


work (1960). 
Parameter Variation Scheme 


The trial wave function can be improved both by selection of 
the best set of linear coefficients through diagonalization of the 
secular equation, and by minimization of the energy with respect to 
the non-linear parameters O. and a. » doe all 2 eae Ng: ile 1eiga slabs 
cussion, the word energy signifies the total electronic energy of the 
desired root of the Secular equation. 

In the earlier modifications of the diatomic program, a pure 
network search was used, because there was no way of knowing at that 
time whether more sophisticated searching routines would get trapped 


id tase maniima: 
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An analysis of the energy dependence on the non-linear para- 
meters has now become possible since large numbers of calculations on 
many types of systems have been performed using the stepwise search. 

Several representative calculations have been analyzed, 
including ground states, excited states, heteronuclear and homo- 
nuclear molecules, atoms, and resonant states. Some definite features 
were observed, leading to the following conclusions: 

1. In true bound states of atoms or molecules, both 

ground and excited, the dependence of the energy 
on the non-linear parameters is smooth. No false 
or pSeudo-minima were found in any of the extensive 
number of calculations checked. 

2. The energy is considerably more sensitive to the 

choice of the 5 parameter for a given orbital than 
the choice of the Qa of the same orbital. 

a. The parwameters*ofr thesorbivals, if originally 
chosen in a "reasonable" fashion (not such that 
the maximum density of a o2s lies inside the 
maximum density of a ols, for example), are 
relatively independent of each other. 

4, Resonant states exhibited false minima in the 8 
parameter of the added electron (outermost orbital), 
and the @ dependence was nearly flat. The energy 
is particularly sensitive to very small variations 


of § in the region nearest the minimum. 
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Thus it 1s now possible to use a more sophisticated parameter 
searching routine, as long as some option is retained to avoid the 
false minima in resonances (or other false minima if they are sus- 
pected). 

Consider the case of a set of parameters chosen so that the 
energy is near the minimum value attainable with the trial function 
being used. Let Xs represent the non-linear parameters, either Oo. or 
Oe BE. =the energy calculated with the set Xs and va) the set of 


au 


increments (x, = a). The energy can be expanded as 
e 
CE iL O° E 
= =: a oS 
B= &, (2 Wi +3 Sx,ox, J ¥1¥j 7 +++ (62) 
: il : sees 
i O al j O : 


Pmaigetne Minimum will be attained if 
OE 
me mae (62) 
1 


Dropping all terms of order greater than two in the expansion and 


substituting the minimization condition, we obtain the set of equations 


OF ae) 
+ ) ls =e (63) 
(5:"] ; ( — 1 


which can be solved for the revised set of parameters if 1/2(n+2)(n+1) 
calculations of the energy are done (n is the number of parameters in 
the set x,). This can be seen to be the equivalent of assuming the 

energy dependence on each parameter to be a parabola, and finding the 


minimum of the set of parabolae which are linked by the cross terms 





EEO 


x;x, (4 = j). Ransil (1960) applied a similar procedure to diatomic 
molecule calculations and found that if the parameters were nearly 
independent, this method converged in three to five passes through 
the set of 1/2(n+2){nt1) energy calculations, using the new y, to 
determine a revised Eo: 

However, if the parameters are independent, or nearly so, the 
effect of the cross terms should be negligible, and an equivalent 
result should be obtained if each of the parameters x, were fitted to 
a parabola holding the other parameters fixed, and the set x, Stepped 
through one by one. This requires only 3n energy calculations per 
pass through the set X,; and as there are two non-linear parameters 
per orbital, produces a nice reduction in the amount of work necessary 
to minimize the energy with respect to the non-linear parameters if 
there are two or more orbitals. 

This process, with a slight modification to be described below, 
is the one used in the present diatomic program. It also converges 
quite rapidly, with two passes normally minimizing the energy to five 
or six figures. 

The actual scheme works slightly differently for variation of Bs 
and Of O,. Since the initial set of & parameters may be quite far 
from the optimum set, each parameter is stepped with a large step size 
in the forward or backward direction until there exists an me Ge) 
which is less than both E. , (6, 4) and Ee oes where i indexes 
the number of the calculation. If the energy continues to improve 


in a certain direction, the step size is increased to hasten finding 
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the 5 4 which causes the energy to turn back upwards. The optimum 


1 
parameter 6 is then determined by fitting to the parabola passing 


through the lowest three points E. ; E and E. In this pro- 


al la als 
cess, more than three energy calculations may be required before 
fitting, but it has the advantage that the first pass through the & 
set always produces a revised BE, which is quite close to the true 
meoamum of the functional. 

since the Q@ dependence is more flat, the search procedure for 
the first pass through the @ parameters only makes three calculations 
per Q, at the values QO,» 1/2 Qs» and 3/2 a: On subsequent passes 
through the @ set, the method is the same as for 8. 

After one complete pass through all the 8 and Q@ which are to be 
optimized, the step size is reduced to one-third the starting value 
and the process 1s repeated. 

In the case of resonances, a Stepwise search over the parameters 
of the outermost electron is still necessary in order to avoid 
becoming trapped in a pseudo-minimum. This option is included in the 
program by setting a tag which both demands a stepwise search in the 
forward direction, and determines the number of steps to be taken. As 
in the previous modifications of the program, the programmer-user 
must apply his experience and physical intuition to determine the 
magnitude of the step size and the number of steps which will locate 
the parameter near the true minimum without becoming trapped. After 


each of the designated outer orbital parameters has been linearly 


searched by this method, the optimum set of stepwiSe parameters is 
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re-minimized by the parabola search method. 

In all of the above procedures, the search may be terminated 
externally by setting of a sense switch. In the parabolic search, 
variation of a particular parameter is eliminated in subsequent 
passes if the program determines that the energy is relatively 
independent of it. This can also be done externally by proper 
preparation of the input cards. 

If the program is not terminated externally, it proceeds 
automatically through the list of parameters to be varied, reducing 
the step size after each pass until it is lower than a preset thres- 
hold step size. At this time, a final calculation is made with the 
optimum set of parameters and whe search for the problem is terminated. 

iInNvciis modification of the program, iG 35 also pocci pseu 
treat the internuclear distance R as a parameter, in that after 
completion of a parameter search at one R value, R can be incremented 
or decremented and the problem restarted. In this case, the non-linear 
parameters are scaled by Beaty eae and the @ search assumes that it 
has already made one pass. 

Thus, if computer time is available, it is possible to calculate 
the entire potential curve for a particular molecule in one run, 
stepping through R and minimizing the energy in terms of the non- 


linear parameters at each desired R value. 
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APPENDIX II. DETAILS OF THE DIATOMIC MOLECULE PROGRAM 


AND INSTRUCTIONS FOR ITS USE 


Overall Logic 


General Comments 


The diatomic molecule program as written is designed for a 
32,768 word high-speed computer with at least four peripheral storage 
units (tape or disc). It is in the FORTRAN IV language, with no 
machine language programs, in order that maximum machine-to-machine 
compatibility be achieved. To the same end, no advantage has been 
taken of special features of some machines, such as reading tapes back- 
ward, multiple entries to subroutines, availability of additional 
peripheral equipment, etc. Variable dimension statements and logical 
IF statements have been used. 

A rather extensive list of options has been included, in order to 
do the most general diatomic molecule problem without changing the pro- 
gram internally, and still allow some special treatment of specific 
problems. 

The program is designed to calculate the energy and electronic 
wavefunction of heteronuclear or homonuclear diatomics. In the homo- 


nuclear case, gerade or ungerade symmetry can be achieved either 


through the use of symmetry adapted one-electron orbitals, or by an 
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option of reflecting the entire molecular function through the 
inversion center and taking the proper linear combination of the 
original and inverted functions. 

The spin portions are set to handle the complete open-shell 
(different orbitals for different spins) problem automatically for up 
to a 28 x 28 set of representation matrices, which includes the full. 
eight electron case. By a change in the input deck, up to twenty 
electrons may be handled with partially or completely closed shells. 
The spin routines are also arranged so that advantage can be taken of 
the reduction in numbers of permutations which can be introduced 
through spatial orthogonalities, by never generating those permutations 
which would make no contribution to the secular equation. 

The one-electron spatial functions are two-centered elliptical 
orbitals. A transformation routine is included which allows linear 
combinations of these functions prior to diagonalization, so that both 
fixed and variable linear coefficients are possible. Space is 
available for up to sixteen functions and a 16 x 16 transformation 
matrix. 

Non-linear parameter search options include a combination net- 
work-gradient search and a pure stepwise network search, with the 
ability to vary one function independently, or two functions simul- 
taneously. Through appropriate tagging, one, both, or neither of the 
two non-linear parameters of a particular function may be varied, and 
the search routine has the ability to set these tags internally if the 
energy is found to be relatively independent of the parameter being 


varied. 
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Through the use of sense switches, the programmer may also 
externally terminate the search, skipping immediately to the final out- 
put routine with the best set of parameters achieved up to that time. 

The internuclear distance may also be varied, with a new calcu- 
lation of the electronic wave function, including non-linear parameter 
Variation at each RK value. 

Thus the program presented here is a compromise among generality, 
speed, efficiency, and size. For selected specific types of problen, 
some of the options may be removed if desired for further optimization 


of these factors. 
Block Structure 


The logic flow of the program is divided into five main branches, 
which are handled by an input and controlling main program. The 
function of each of these branches is as follows: 
Branch l - SETUP Branch 
Reads auxilliary tape identifiers, loads permanent data 
tables for integral calculations, makes spin-permutation 
tape. Re-entered only if spin problem is changed. All 
output from this branch is labeled COMMON blocks or on 
tape (disc). 

Branch 2 - INTEGRAL Branch 
Controls calculation of all one- and two-electron spatial 
integrals. Makes integral table tape of indices, parameters, 


one-electron integrals, and G, ,{x,p). Combines all G, 
Lk oA 
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with Git to make tape of all two-electron integrals. Upon 
exit, all one- and two-electron integrals are on tape for 
Subsequent iterations, and in fast memory for combining with 
Spin. 

Branch 3 - MATRIX Branch 
Applies transformation to spatial integrals if required. 
Combines spatial integrals with spin matrices from tape to 
make total secular equation. In homonuclear case with total 
molecular inversion option, secular equation is block 
diagonalized in this branch. 

Branch 4 - SEARCH Branch 
Solves secular equation. Varies non-linear parameters as 
desired, prints intermediate output. If parameter search 
not completed, return to INTEGRAL branch, otherwise proceed 
to OUTPUT branch. 

Branch 5 - OUTPUT Branch 
Prints final output, including total secular equation, 
eigenvalues, eigenvectors, converted energies and spatial 
configurations. If program was cut externally by sense 


switch, punches new input deck for convenient restart. 


Concerol of Storage and Internal Overlays 


Of prime importance is the control of the working space inside 
fast memory. In the various branches, and inside the branches, data 
is overlaid when it is no longer needed. In order to prevent des- 


troying information which will be required later in the calculation, 
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an understanding of the places and times that these overlays occur is 
mandatory if modifications are to be introduced, if it is desired to 
use certain routines as "black boxes" for entirely different calcu- 
lations, or if the program is to be updated. 

To facilitate bookkeeping, all fast storage is in four labeled 
COMMON blocks. A description of the contents of these blocks and the 
overlays which occur in them follows. 

/SPACE/ - 309 words. Contains most problem definition data as 

described in input deck summary below, plus 

ITX - iteration counter. l-first pass; C-intermediate passes; 
3-final pass 
FPARD - best set of non-linear parameters, stored from FPAR 
as determined 
EO - best energy, corresponding to calculation with FPARG 
NO - tag for parameter being searched. 1-delta; 2-alpha 
MO - function number of parameter being searched, i.e., 
search over FPAR(NO,MO) 
ICTAG - 64 locations available for additional tags as desired. 

Data in this block is generally permanent, with the only overlays 

being resetting of tags, changing parameters, etc. 

/DATBLS/ - 759 words. Permanent information. Not changed during 

run. See below for details. 

/ Otay - 268 words. Transformation matrix and tags. Matrix is 


recalculated if NOTRF > $, but not overlaid by other data. 
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/WORK/ - 12,518 words. Main working area. Overlaid by all 


branches as follows: 


i 


SETUP - overlays controlled by subroutine SPIN. Used to 
contain all information for generation of spin tape. 
Completely destroyed upon exit from spin routine. 

INTEGRAL - during making of one-electron tapes and com- 
bination to make two-electron tape, /WORK/ is composed of 
two arrays JTA and JTB, dimensioned 146@x4, and the one- 
electron integrals Sl and Hl, singly dimensioned at 

136 words each. Sl and Hl are stored as calculated. The 
two-electron tape is read upon completion, storing the H2 
integrals over JTA and JTB. rc JTB is used as temporary 
storage by the integral routines ONEL, BANDI, and GLIST. 
289 locations available in array named DUMBO. 

MATRIX - at entry, contains H2 (9317), SPRI, HPRI (Sl, Hl 
of INTEGRAL branch), transformed one-electron sets renamed 
S1,H1l, one 28x28 spin matrix, one permutation, and room for 
the total S,H stored triangularly (1806 words). At exit, 
H2 is overlaid by the total secular equation OVP, HAM 
which is S,H unpacked and stored as full two dimensioned 
arrays. Array named DOM has 19 locations available. 

SEARCH - six 4exke matrices for secular equation and eigen- 
vectors, including two copies of secular equation, one for 
printing if desired, and one which is destroyed during 


diagonalization. Several vectors dimensioned 42 for 
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eigenvalues, etc. Array named FINK has 1760 locations 
available. 
5. OUTPUT - essentially same as SEARCH, except FINK contains 


only 1672 locations. 


famitations 


The following limitations must be observed to prevent premature 
termination of the calculation by built-in error checks: 
Maximum number of one-electron functions 16 
Maximum number of electrons in IPER jad ©, 
Maximum number of spin eigenfunctions @ 28 


Maximum number of spin functions 2 -- 


Maximum dimension of secular equation he 
largest n,, m, in 
Largest Qa, 19.9 
smallest 5. oe 


Error Messages 


Error messages may occur during calculation if the internal 
checks are positive. Some of these errors will not cause premature 
termination of the program, and others will. The following is a sum- 
mary of these error messages, their causes, and the action taken by 
the program: 

1. WRONG PERM TAPE xxxxxx XXXxXXxX 

The identification tag on the spin-permutation tape does not 


agree with the tag required by the formula 
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JCODE = 1O00*N + 10*10SH + MULT 
calculated when the problem is defined. This can be due to 
a mis-positioned tape or to having mounted the wrong tape. 
The program repositions the tape and either continues or 
exits, depending on the condition of a recheck. The first 
number in the message is the desired number, and the second 
is the number read from the tape. 
SEC EQN TOO BIG xxxx 
The dimension of the secular equation requested by the input 
parameters (NCONXND) is larger than 42. 
PNDICES VON CHECK xX XX XX X°5xXxX xx 
Subroutine INTEG has detected that the indices on tape KR 
and the indices in the integral calculation loop do not 
correspond. The indices Ll, I and J are printed in pairs 
with the supposed Ll, I, and J read from tape KR. 
BANDI - NO CONVERGE 
The uw summation of Eq. (60) did not converge in 36 terms 
or less. 
BANDI - PARAMS OUT xx xx xXx xx 
One or more parameters are out of range. The four numbers 
printed are T, m, lol, and M, where T, m, and | o| are defined 
in Appendix I and M = 28 + /al/1.95, which defines a dimension 
and cannot exceed 60. 
GLIST - GENERAL ERR 


self-explanatory 





LO; 


del. 


pled 


GLIST - PARAMS OUT xxxxxx xx xx 

Parameters out of range in subroutine GLIST. The three 
numbers are 5, n, and IM = I4 +n + Ot + 2, where 8 and n 
are explained in Appendix I, and IM is a calculated 
dimension which cannot exceed 45, 

ENERGY. LT. MIN 

The calculated energy is lower than the experimental or test 
energy punched on the title card. An exit is made to the 
OUTPUT branch, after which a new calculation may be started 
if more input cards are present. 

GR.SRCH SEEKS MAX 

The gradient or parabola search is seeking a maximum rather 
than a minimum. This is an information message only. 

END SRCH OF FCN xx xx 

The energy has been found to be independent of a particular 
parameter, and the tag has been set to discontinue the 
search of this parameter. Its indices are printed for 
information. 

OVRLAP NOT POS DEF xxxxxxx xXx xx 

The overlap ates is not positive definite, If it is 
positive semi-definite, this is an information message only. 
If the magnitude of the negative diagonal element is less 
than 10° |, the program assumes that the matrix is numerically 
positive semi-definite, deletes the row and column cor- 


responding to this element in both the S and H matrices, and 





A short summary of the programs follows. 


continues. 
easily be seen. 


than 10. 


7 


This results in an eigenvalue of 107° 


dee 


which can 


If the negative diagonal element is greater 


in magnitude, the program makes an error exit. 


Specific Programs and Subroutines 


Those routines which ° 


are marked by asterisks are discussed in detail below the summary. 


The logic of the unmarked routines can be followed without difficulty 


in the program listings themselves. 


HUY Old 
DIAT 
SETUP 


TABLE* 


SEIN* 


CONFIG* 


HDG 


INTEG* 


ONE L* 


BANDI 


GLIST 


TWOEL 


COMB IN 


A Gadd 


M 


BRANCH 


MAIN 
SETUP 


SETUP 


SETUP 
SETUP 
MAIN/ 
OUTPUT 
INTEGRAL 


INTEGRAL 


INTEGRAL 


INTEGRAL 


INTEGRAL 


INTEGRAL 


PURPOSE 


DATA INPUT & MAIN CONTROLLING PROGRAM 
INITIALIZATION OF TABLES, ETC. 


FIXES PERMANENT DATA BLOCK, 
athe 


POINTS, WEIGHTS, 


PREPARES TAPE OF SPIN- PERMUTATION MATRICES 


READS CONFIGURATION TABLES INTO LCON (config. 
array ) 


PRINTS TITLE PAGES IN OUTPUT 


CONTROLS PRODUCTION OF ALL SPATIAL INTEGRALS 


MAKES ONE-ELEC TRON INTEGRALS FROM BANDI 
OULEUL 


MAKES B,A,AND I INTEGRALS 
MAKES GA 
INTEGRALS 


x) FROM & INTEGRATION AND I 


7 WITH Gay 


CONTROLS CALLING OF TWOEL FUNCTION 


COMBINES G, TO MAKE Cllr ail) 





Peck LLP 


TAPE 


MATRIX* 


TXFORM 


TRANS 


TRX 


SEARCH 


D23CHO 


KIG 


OUTPUT 


S 


S 


S 


S 


BRANCH 
INTEGRAL 


MATRIX 


MATRIX 


MATRIX 


MATRIX 


SEARCH 


SEARCH 


SEARCH 
OUTPUr 


OUTPUT 
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PURPOSE 
READS AND WRITES INTEGRAL TAPE 


COMBINES SPACE AND SPIN PARTS TO MAKE TOTAL 
SECULAR HQUATION 


PREPARES TRANSFORMATION MATRIX FOR SPATIAL 
FUNCTIONS 


APPLIES TRANSFORMATION TO ONE-ELECTRON 
INTEGRALS 


APPLIES TRANSFORMATION TO TWO-ELECTRON 
INTEGRALS 


ENERGY CALCULATION AND PARAMETER SEARCH LINK 


SOLVES GENERAL EIGENVALUE PROBLEM AND ORDERS 
ROOTS, VECTORS 


JACOBI MATRIX DIAGONALIZATION 
PRINTS FINAL OUTPUT 


PRINTS MATRICES AS SQUARE OR TRIANGULAR ARRAYS 


*TABLE - FIXES PERMANENT DATA BLOCK, POINTS, WEIGHTS, ETC. Sub- 


routine TABLE generates 


definitions are: 


and 


Die 
ae) 


DDE. 4 
1, J 


DDe. . 
yd 


Dl, DDL, and DDe2 for BANDI and GLIST. Their 


= (i-1)1(2j-1)/[(i-g) 1! (atg-1) 19) (64) 
= [(itj-2)1(2i-1)/(4-3) 1]? (65) 
= ((itj-1)(4-g41)/(44°-1)] , (66) 


and they are calculated recursively. 


The numerical quadrature points and weights are read in this 


routine. 


They are read as data rather than loaded as a permanent 
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BLOCK DATA Subroutine to allow the use of several different quadrature 
formulae, with different numbers of points. Space is available for up 
to a forty point integration raule. The integral limits (l,) are 
transformed to (0,1), so the points are in increments of 1/x and the 


weights are 
il 1 
wer an 
Ww, = x,w, 2/ (x, - 1)? (67) 
where 


Le 


x, is the 4h oint in POINT 


is the value punched on the input cards, WT(K) 


Ws is the tabulated weight from the numerical quadrature rule. 

Card decks have been prepared and are available for the original 
34-point Newton-Cotes formula used by Harris and Taylor (1963, etc.), 
and for 16,20,24, and 32 point Gaussian formulae. Table 9 gives a 
comparison of two-electron integrals for various 8 values using these 
formulae. The 34-point formula should be used in any check calcula- 
tions which attempt to duplicate the results of the old Taylor-Harris 
or Taylor calculations. 

*SPIN - PREPARES TAPE OF SPIN-PERMUTATION MATRICES, Subroutine 
SPIN applies the methods described in Appendix I, (84 ff.) to generate 
the permutations IPER, and their NDXND representation matrices 
UMAT(I,J). <A given permutation (read or calculated) is applied to 
the spin product functions ISP, and an ordering vector ISSPM is created 
such that if the permutation is denoted by P, and P on ISP(I) yields 


ISP(K), then ISSPM(K) = I. Using this ordering vector, 
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UMAT(I,L) = - USET(I,J) + USET(L,K)|-SIG 
J=1 
where USET is the spin eigenvector matrix, K = ISSPM(J), and SIG = +1, 
the parity of the permutation P. 

*CONFIG - READS CONFIGURATION TABLES INTO LCON. Subroutine 
CONFIG is simply a reading routine using variable dimension statements 
which allows maximum use of the 150 words allocated for the configu- 
ration table LCON. LCON is dimensioned as an NCONXN matrix rather than 
with fixed dimensions, so that (5 two-electron, 50 three-electron, etc. 
configurations are possible. CONFIG reads, prints, or punches LCON 
in this type of storage. Note, however, that the maximum dimension 
of the total secular equation is 42x42, and that NCONXND may not exceed 
fais Limit. 

*INTEG - CONTROLS PRODUCTION OF ALL SPATIAL INTEGRALS. The 
second major branch is controlled by subroutine INTEG, which handles 
the calculation of all spatial integrals. This is done in two steps, 
the first being the calculation of the one-electron integrals Teng 
ee and oe forsale aad: 4) (iS j),» and second, taking all possible 


combinations of a with G., to make the two-electron integrals. 


a 
To save storage space, and due to the symmetry of the spatial 
one-electron matrix elements, the oe and aera sets are stored as 


vectors containing only the len(nt1) elements of the lower triangle. 


The (i,j) element of the matrix can be looked up by 


INDX(I,J) = MIN(I,J) + [(MAX(I,J)-(MAX(I,J)-1) ]/e2 


where MIN(I,J) is the lesser of i and j, and MAX(I,J) the greater. 





ey. 


By the same token, a double application of this method allows 
storage of the two-electron integrals in a single vector H2, containing 
1/2 m (m+1) elements where m=1/2 n (ntl). In this case, the > inte- 
grations are done later. 

A single (i,j) integral table requires 1460 words — 19 words for 
indices, parameters, and one-electron integrals, and the (40x36) 

G5 004) matrix (also treated as a one-dimensional vector by INTEG — 
unpacked only in GLIST and TWOEL). 

For optimization of input-output time, since these tables (named 
JTA, JTB, or ZGl in the routine) must be stored on tape or disc, they 
are generated in blocks of four, and written on tape KW each time four 
sets have been completed. One to three blank tables may be written 
in the last block of four in order to make the total number of tables 
@ multiple of four. If a previous integral tape KR exists (i.e., 
second and subsequent passes of parameter variation scheme), the 
pertinent blocks of tables are read from this tape before actual 
integral calculation begins. 

Subroutine ONEL is called for generation of each table. Upon 
return, the ay and oe are extracted, normalized, and stored in 
fast memory for subsequent branches. 

The second, combination, stage of INTEG uses COMBIN and TWOEL. 
Every block of four tables must be combined with every other block, 
making ten or sixteen two-electron integrals depending upon whether 


a block is combined with itself or another block. 
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To effect this combination of Gok tables from the tape, the 
following logic is used to minimize tape-handling time, which 
experience has shown to be the rate-controlling step in the entire 
energy calculation. 

Tet there be NTOT blocks of four tables. The last (wrot’*) 
block is still intact in fast memory. This block is combined with 
itself while the tape is rewound. The tape is then read forward and 
blocks 1,2,3....., NTOT-1, are combined with NTOT. At this point, 
the wrort’? block has been completely combined, so NTOT is reduced by 
one. The tape is now read backwards (actually, two BACKSPACEs followed 
by a read to remain in the FORTRAN IV language), combining the new 
wror’? block with all the blocks down to block one. This process is 
continued back and forth dropping a block each time as it has been 
fully combined until all the blocks have been eliminated. For example, 
if there are five blocks of four tables, the logic would produce the 
following order of block combination: 

55s oO, SS NOM yy 3 ee leone Gres een 

Since the single index of each two-electron integral is calcu- 
lated from the i, j,k,£ indices in the tables and written on the two- 
electron integral tape KINT along with the integral, this peculiar 
order of calculation has no effect. 

Upon completion of the combinations, the two-electron integrals 
are read into fast memory for the next branch. 


*ONEL - MAKES ONE-ELECTRON INTEGRALS FROM BANDI OUTPUT. Sub- 


routine ONEL is called by INTEG with a particular i and j, and a 


designated table of 1460 words. ONEL compares the parameters in the 
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table (if an old table exists) with the current set and either returns 
if there has been no change, or calculates a new set of one-e_ectron 
integrals and G vector through calls to BANDI and GLIST. 

The layout of each individual table is as follows {Greek syrbols 
as defined in Appendix I ( ff.)]: 
WORD DEFINITION 


ik NLP(1) - v.-v. = a 


2 niP(2) - (lv, +lvi-fol)/2 = + 
3 NLP(3) - m, +m. = m 
4 NLP(4) - number of steps necessary for convergence or I integrals 
(i.e., upper limit of » summation’ 
5 NLP(5) - nj tn, 35a 
6 NLP(6) - m.-m, = m/ 
1 og 

= a = rg 

a NLP( 7) n; n, n 


8 NOTA(1) - I 
9 NOTA(2) - J 
10 EPOT - potential energy term, independent of & 


11 INDX - J+I(I-1)/2 (J < I) 


12 Maximum length of table - 19+JMAX*(NLPM4) - o 3 
13. FLP(1) - 5, +6, = 6 

= = = / 
Lie ChiGPee)) oF 8 a) 
15> FLP(3) = ata, =a 
16 FLP(4) - a.-a, = a/ 

ed 

Li aee St lo) 


18  EKIN - kinetic energy term, independent of R 
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WORD DEFINITION 


19 GINY - point at infinity for Giy zero except in 34-point Newton- 
Cotes integration (see Harris, 1960 - Eq. 115) 
20 Gs - integral table, of length JMAX-(NLP(4) - |o|) 
*MATRIX - COMBINES SPACE AND SPIN PARTS TO MAKE TOTAL SECULAR 
EQUATION. Subroutine MATRIX calls TXFORM and TRANS to prepare the 
Spatial one-electron integrals in final form. 


The spatial integrals and the spin matrices are combined to 


make the secular equation 


> U(P) [H(P) - ES(P)] = 0 
je 


storing H(P) and S(P) as row vectors of dimension 1/2 NMAT-(NMAT+1), 
where NMAT is NCONXND. 

the routine makes one pass through the spin tape, preparing all 
contributions of a particular permutation matrix to each matrix ele- 
ment in H and S before reading the next set on IPER amd UMAT from the 
tape. 

The function TRX is used to look up the proper two-electron 
integral and apply any transformations to it, plus doing the p inte- 
eration. It is called withthe indices 1,j),k,2"and retUmms vie imte- 
gral ist 1) sees |eg) e 

After the secular equation is complete, it 1s unpacked and stored 
as two full matrices, HAM and OVP, for use by the diagonalization and 
search routines. 

In the homonuclear case, if the full molecular inversion option 


has been selected, the secular equation 1s pretransformed at this 
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point into block diagonal form, with two blocks of equal size. The 
block corresponding to the desired symmetry (g or u) will be in the 
upper left, and only this block will be carried over into the next 


branch. 


Input Deck Setup 


The input deck is divided into three sections: 

1. Points and weights for numerical quadrature (SETUP branch) 

2. Spin information (SETUP branch) 

3. Problem definition and electron information (MAIN branch) 
Section 1 is common to all calculations in which a specific quadrature 
rule is desired, and is retained in memory intact throughout the entire 
run. Section 2 is common to all calculations involving the same number 
of electrons and a specific spin multiplet. It is used to produce the 
spin tape, which is retained until subsequent data cards redefine the 
spin problem. Section 3 is only common to the calculation which is 
immediately in progress. 

One comment is necessary here. This program was designed to be 
compatible with any machine which accepts a FORTRAN IV language. Due 
to differences from system to system, all input-output statements 
involving tapes (common input tape, work tapes, etc.) use symbolic 
names. The actual numbers can be read from data cards, with one 
important exception. The symbolic tag of the common input tape for a 
particular system must be set internally. Therefore, even though this 


is not part of the input deck, this is included here. [In the main pro- 


gram, the variable IN must be set to correspond to the logical tape 
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number for the system common input tape. (Example - for the Honeywell 
Automath 800/1800, IN = 2, for IBM 7094/7044 pcos, IN = 5). 


The setup for the input deck proper is as follows: 


BRANCH CARD NO. FORMAT Cit Seh 
SETUP cl lel KW, KR, NT,KINT, IOUT, NPCH 


KW - Write tape for integral tables 

KR - Read tape for integral tables 

NT - Spin tape 

KINT - Two electron integral tape 

IOUT - Logical number of system common output tape 


NPCH - Logical number of system common punch tape 


SETUP 2 ( 6X, 2F8, E30. 12, 2A6) JMAX, NMINUS, WTINF, RULE 


JMAX - Number of points in quadrature formula 


NMINUS - Number of points to be accumulated negatively 
(See WT below) 


WIINF - Weight of point at infinity (= 0 for Gaussian quadratures) 


RULE - Alphabetic name of quadrature formula (limited to 
12 letters) 


SETUP 3-12 (6X, 4H18. 12) POINT 
POINT - Array of quadrature points (Maximum = 40 points) 
SETUP 13-22 (6X, 4E18.12) WT 


WT - Array of weights corresponding to points. These contain 
the square root of the weight given by the quadrature 
formula. Since some may be negative leading to imaginary 
roots, only absolute values are stored, with the first 
MINUS weights being the negative ones. 


oBTUP pe (1213) JNT, NEL, NISP, MULT, ND, I[OSH, JGEN 





JNT - Tag. JNT=O, make new spin tape; JNT=1, spin tape already 
mounted 


NEL - Number of electrons 


NISP - Number of spin product functions (eigenvalues of Ss.) 
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BRANCH CARD NO. FORMAT LIST 


MULT - Spin multiplicity 2S+1 


ND - Number of eigenfunctions of go (= dimension of a one- 
configuration block of the total secular equation) 


TOSH - Number of open-shell electrons (See App. I and II) 
JGEN - Number of hand-generated permutations (See App. I and IT) 


SETUP 2(deck) (1415) HSE 


ISP - Spin product functions. Read as column vector of 
dimension NISP 


S20UP 3( deck ) (5F14.10) USET 


USET - Coefficient matrix for eigenfunctions of a Dimension 
NDXNISP, read as column vectors. 


SETUP 4 (deck) (F9.4, 2113) SIG, JPRI, KPRI 


SIG - Parity of a hand-generated permutation 





JPRI - Hand-generated permutation (Example - +1. 3 4 1 2) 
KPRI - Aux. permutation with double occupancies noted 


Note - if JNT = 1, card decks 2, 3, and 4 are omitted. 


MAIN ak (8A6,3F10.6) TITLE, Al, Ae, EFIN 


TITLE - Run identification of 48 characters, to be printed on 
each output page 


Al - Atomic energy of atom l 


A2 - Atomic energy of atom 2 (Al,A2 used to calculate dissoci- 
ation energy) 
EFIN - Experimental energy. If E ale < Bean 
assumes an error exists and stops. If no error test is 
desired, set EFIN = oO. 


MAIN a (7F10.4) ZA, ZB, RMAX, STEP3, STEP2, RMIN, RSTEP 


ZA - Nuclear charge of atom A 


» the program 


ZB - Nuclear charge of atom B 
RMAX - Maximum internuclear separation. 
STEP3 - Maximum step size in variation of non-linear parameters 


STEP2 -~ Minimum step size in variation of non-linear parameters 
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BRANCH CARD NO. FORMAT LARS b 


RMIN - Minimum internuclear separation. If blank, only RMAX 
is used. 


RSTEP - Decrement to RMAX until R < RMIN 


3 (1273) N,NHUG, NCON, NFCN, NPC, NE, LSRCH, 


NOTRE, NU, KOLD, IS 


N - Number of electrons (used as check on NEL and spin tape) 


* 


NHUG - Tag. O = Heternuclear; +1 = Homonuclear with symmetrized 
orbitals; te = Homonuclear with inversion of total 
molecule; > 0 = gerade, < 0 = ungerade. 


NCON - Number of spatial configurations 
NFCN - Number of spatial functions 


NPC - Number of functions whose parameters are to be varied. 
vor single caleulation,, set 0. 


NE - Bigenvalue to be minimized. 0O = ground state, l = fee 


Crereed Stave, etc. 


LSRCH - Search tag. 0 = standard step and gradient search. 
> © = linear search of LSRCH steps in positive dineerier 


NOTRE =- Spatial transformation tag. < 0, read C.. from cards; 
= Of no transiormat ion; = 1 projyece out’ first NU 
funetions from specified orbitals; = 2, orthonormal ze 
all spatial functions. 


Hi== see NOTRE. If NOTRE = 0, NU can be used to start parameter 
search in NU + 1 spatial function. 


KOLD - Variation tag. = 0, vary each spatial function inde- 
pendently; = 1, vary i and i + NFCN/2 simultaneously. 


ic = Search tag. 0 = search O as 3/2°O, Oy ecm eo ccomene 
as 5. (Special note: if an old integral tape KR exists, 
iS] 2 will set [fx = 2. This provides a cemventent 
restart if the last KW tape has been saved.) 


4 (deck) (15, 2F10.6,613) I,FPAR(1,1),FPAR(2, 1), (NFPAR(J, 2), 
J=1, 6) 


Spatial function cards 
IV Minetion- number 
BEAR Gly) aes 
FPAR(2,1) - a, 

7G bee Beas ry 





BRANCH CARD NO. 


NFPAR(2, 1) 
NF PAR (3, I) 
NF PAR (4, T) 


NFPAR(5, I) 


NFPAR(6, I) 


5 (deck) 


USS. 


FORMAT EESG 


Ta 
ae 


Vv. 

1 

Parity Of symmetrized Orbital 0 —"serade.] 15 — 
ungerade 


Projegpion tan. —"O = NO preyeculon, — lb. pro)jccr 
Ollie ik: NU functions 
Variavion tag. — Oy "search Dou Gvdnd 75 
search Q only; = 2, search 5 only: | — espe eccomcn 
or This-funecuien 

G2 I, (LCON(1I, J), J=1,N) 


ie scOurieuracvion dumber 


HeON = Array of functionmmumbers belonging Vembne conristnmavvods 


6 (deck) 


Geo ((c(1,7), J=1,N), 1=1,N) 


C - Transformation matrix for spatial functions. No deck unless 


NOTRF = 


~1 
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APPENDIX III. ILLUSTRATIVE EXAMPLES 


AND HINTS - THE NO* CALCULATION 


Illustrative examples are furnished to allow the reader to see 
just how an input deck is set up for a particular type of problem, 
and as calculations which can be used to check that the program is 
running properly on a particular computer after the minor modifi- 
cations necessary to make it compatible with that computer have been 
made. 

Three problems are listed in the sample input decks of Table 
The first of these is a LiH calculation, illustrating the hetero- 
nuclear capability of the program. This is a one-configuration four 


electron singlet calculation, and the correct results for the lowest 


meet are: 
1. Total electronic energy = -9.708228 au 
2. Electronic + nuclear (potential curve) = -7. 708228 au 
3. Dissociation energy = +7. 708228 au 


The input deck for this run is self-explanatory once the summary 
in Appendix II is understood. 


The second example is an H, calculation using the symmetrized one 


2 


eleetron orbitals and two configurations, C-c_iso_is’ + Cro liso lc 
ee g ou u 


This is for the ground state and the calculated energies should be: 
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-1.869461 au 


(ieee letar oLectronic enersy 


2, Electronic + nuclear (potential curve) = -1.155176 au 


3. Dissociation energy = Jip iaie wae 
Note that NHUG = +l for a gerade state with symmetrized orbitals. 
thard, the Hy calculation is repeated using the option of total 
jmeaersion Of the molecule through the inversion center, selecting the 
gerade linear combination ¥ + i¥. Several things should be noted 
here. The first is that the reflected counterpart of each function 
(same function with a replaced by -a) is included explicitly in the 
immts deck, and that all unreflected functions precede Une rei lecred 
ones so that if there are m unreflected functions 1,2,...m, and each 
Senresponding reflected function is denoted Me ener the input 
MubetIOns are numbered 1,2,....,m, mtl = sie idle peta anne Oconee » mim=m. 
Hae contiguration cards are treated similarly, in that cach comia2r— 
RaclOumcards are treated samilarky, in@that eCdehmecmirsuralion madcue. 
unreflected orbitals must be paired with a configuration composed of 
tiescorresponding retléected orbitals, wath all unrerlectedseen— 
Piearaelons. lie program combines the reflected Contieural once wae 
rotation by nm and adds or subtracts the®paired cenfisurations as 
necessary for g or u. Note that NFCN and NCON are set for the total 
number of functions and configurations. NHUG = +2 is the tag which 


instructs the program to perform the rotation and block diagonalize 


the total secular equation taking only the g block. 
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ihemenecrpies in this case are; 


1. Total electronic energy == OlOoveau 
2. Electronic + nuclear (potential curve) = -1.076781 au 
3. Dissociation energy = -.076781 au 


Note that although the same functions were used in both Hp 
calculations, the one with symmetrized orbitals has a slightly lower. 
energy. This is because the symmetrized orbitals allow two linear 

eee sticients, while in the total molecular inversion option there is 
emivecne. The total inversion, in effect, gives C, (co+0/07) = 

eo ou)? and we see that C) = Vf 2/2. However, the total 
inversion requires only 3/2 to 3/4 the work required by the symmetrized 
orbitals, and if only ee were used in the symmetrized case, a poorer 
answer would result. 

The ieee + Ciel calculation can be improved by varying the 
parameters of the oF orbitals, Samee they are nomemcstricteds ome 
equal to the parameters of the oe orbitals as iS Ghe case of wera 
inversion. They were set equal purely to demonstrate the relative 
effect Of the two options. 

Pome specific hints are also helpiuk for actual use ot sue 
program on particular computers. 

On machines having some type of CHAINing feature, the programs 
themselves can be overlaid. For example, the IBM 7094/7044 DCOS IBSYS 
monitor, version 13, was used with three links. The INTEGRAL, MATRIX, 


and SEARCH branches were included in one link for optimum efficiency, 


and the SETUP and OUTPUT branches formed the other two links. 
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An additional space-speed saving can be effected in systems 
with adjustable input-output buffer sizes, such as the Honeywell 
Automath 800/1800 monitor, by buffering the two-electron integral 
mape tO the smallest possible size, since only two words are written 
mer iogical record; and by buffering the spin-permutation tape to the 
emealest size which will include N + NDXND words in one logical record: 
Tapes KW and KR should be buffered with the largest double buffers 
possible, since tape handling by subroutine INTEG is the time- 
fomprolling factor. 

On the Honeywell 800/1800 system, a small calling program is 
Mececsary to set up the monitor for use with variable Taperidencrt rec. 
The entire deck setup is as follows: 

*JOBID, JOBNAME, 10424345, A8SNGG25, B9DBGG19, GEDBG24G, DTDBG2GG 
[ero jecy number xxx, USe€rs name, any Speciale instruculoncs.o 


* be typed at console for operator. (64 columns/card) 


TITLEMOLEC 
CALL DIAT 
I=0 
WRITE(6)I 
WRITE(7)I 
WRITE(8)I 
WRITE(9)I 
END 
TITLEDIAT --- entire diatomic program deck as presented in App. IV 
*DATA 
6 7 8 9 3 5 (Tape identifiers, corresponding to *JOBID card 
above) 
-Remainder of input deck- 
SEL) 
*ENDFILE 


On the other hand, when the IBM 7094/7044 DCOS system is used, 


the monitor sets up all buffers when variable tape names aré used, and 
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rather than tricking the monitor with the four statements WRITE(n)I 
above, one must use a MAP program which deletes the unused buffers. 
Sepics Ol this program are available on request, or the following 

isomer nay bemused to Create the program for a specific installation. 


To delete the buffer for peripheral unit nu, the MAP program is 


$IBMAP BUF. NODECK 
ENTRY  .UNnn. 
AUN. PLE 
END 


The ENTRY and PZE instructions are repeated with different 
numbers gl, 02, .«. as desired to deleve any wnnecessanveunaine 

ine format for centrol cards, accounting cards, crc uw ecammec 
found by consulting the appropriate IBM manuals and the notices 
published bvueene par prcular ins tallarion, 

im the prozram listings (App. IV) the reader will note several 
comment cards relevant to particular machines, and accompanying 
FORTRAN statements. These must be included or pulled as necessary. 


The comments are considered self-explanatory. 
The NO? Calculation 


As an illustration of the use of the spin-permutation program 
options, the following sample calculation on Not is presented. The 
approximations used were picked more to illustrate the use of the pro- 
gram than for their physical validity and thus the calculation is by 
no means a definitive one. 

To do a complete 14 electron spin problem, 14! 429x429 matrices 


U(P) would be required, which is patently impossible. The approxi- 
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mations used to reduce this number to a manageable one are: 
i. Thewelfosing of the shells whweneincludeseie tén inner 
electrons so that the spatial configuration 1s 


appLoximablely represenved by 


Q a C 2 2 
KP, ~A2P, on2p_ynep__(ols)y (ols), (02s), (02s) 5 (Gi on 


and thus requires only two spin functions 


930, (1/W2(ab-Ba) 1° and 6 vol/W2(a6-pa) }” 7 

2, Orthogonalization of the closed shell (doubly occupied) 
orbitals. This results in a spatial configuration of 
seven shells of two-electron pairs which are pairwise 
orthogonal. 

With the open shell electrons listed first, the auxiliary 
permutation vector for the identity permutation is 

1 2 3.4 55 66 77 88 9a 
and obviously no permutations are necessary between electrms 55, 66, 
etc. which reduces the number of permutations necessary for spin 
reasons to 14! /2?. We now take advantage of the spatial ortho- 
gonalities which were purposely introduced. 

As pointed out in App. I, any permutation which produces three 
or more integrals (i{j) which are zero can be eliminated, since it 
cannot contribute to the secular equation. Inspection of the 
auxiliary permutation vector above shows that only 2O permutations 
among the first four orbitals can contribute to the secular equation, 


and that no more than one interchange of the closed shell functions 


is allowed at a time. 





1he 


Thus, there are 31 allowed interchanges of the closed shell 
orbitals with themselves or with the open shell orbitals 1-h. 

The twenty allowed permutations of the open shell set result in 
620 permutations which can contribute to the secular equation within 
this particular closed shell and orthogonality system, and the result 
is 620 2 x 2 matrices U(P) which is obviously a considerable improves 
mMenueover the complete open shell treatment. 

To make the input deck slightly easier to prepare, the ortho- 
gonalities in the set of four open shell electrons are ignored, and 
all 4! permutations among them are produced in the program. Hence 
there are 31 x 4! = 744 permutations and their corresponding U(P) 
foroduced. | 

Me inpuc deck for this calculation is shown in fable iy wicca 
Pedestie paired orbitals are replaced by their canonical numbers, 
and that IOSH = 4. Since the permutations among the ten closed shell 
Pe-eerons are Explicitly done by hand, no auxiliary permutarien 
Meccoreussneeded an the input deck. The tape adenti7 ierseand 
Wiadracure points have been omitGed. Preparation of Lhewspitm cape 
required about two hours on the Honeywell 800 and about eight minutes 
on an IBM 7094/7044 direct. coupled isystem, 

Onee Tie spam tape was prepared. Tne Not calculation itself 
took about three minutes on the 7094 and about twenty minutes on the 
H-800. The resulting energy was -105.93 au, which is of the right 
order of magnitude, but does not even show binding over the separated 


atom calculations of Clementi (1964). Due to limited time available 
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to the author, no attempt was made to improve this calculation by 
Tearlazdon Of the non-linear parameters or by making more realistic 
physical approximations. Upon completion of the single one-shot run 
dt was decided that the original purpose of this study had been 
Penpeved, that of demonstrating that a calculation on @ large many-— 
electron diatomic molecule is possible and practical with the present 
me@ititcation of the diatomic program. 

The short time required for the actual calculation shows that 
parameter variation and even some configuration mixing are possible 
in reasonable amounts of computer time, and a much improved energy 


obtainable even with the limited spin set used. 
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